


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1996-06 


Analysis of steam and hydronic compartment 
heating systems aboard U.S. Coast Guard 140 
foot WIGB class cutters 


Hurley, James Thomas. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/32089 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 





IDTITC QUALITY tysPtoren 2 


. 1 
at ty 
wht fe 


THESIS 


ANALYSIS OF STEAM AND HYDRONIC 
COMPARTMENT HEATING SYSTEMS ABOARD 
U.S. COAST GUARD 140 FOOT WTGB CLASS 
CUTTERS 


by 
James Thomas Hurley 


June, 1996 


Thesis Advisor: Ashok Gopinath 





Approved for public release; distribution is unlimited. 


WSOTOOT O87 








~ DISCLAINER NOTICE 






THIS DOCUMENT IS BEST 
QUALITY AVAILABLE. THE 
COPY FURNISHED TO DTIC 
CONTAINED A SIGNIFICANT 
NUMBER OF PAGES WHICH DO 
NOT REPRODUCE LEGIBLY. 














REPORT DOCUMENTATION PAGE 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction, searching existing 
data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate 
or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information 
Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork 
Reduction Project (0704-0188) Washington DC 20503. 


June 1996 Master’s Thesis 
SYSTEMS ABOARD U.S. COAST GUARD 140 FOOT WTGB CLASS CUTTERS 
(6. AUTHOR(S) JamesT.Hurley — —“—s—s—s—“—“<‘LSsi‘(CSstsi S—<(Cs—‘CSs—SCSCCCCCS 


PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 


Naval Postgraduate School 
Monterey CA 93943-5000 









8. PERFORMING 
ORGANIZATION 


ig 






REPORT NUMBER 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 








11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect 

the official policy or position of the Department of Defense or the U.S. Government. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited. 

13. ABSTRACT (maximum 200 words) The compartment heating system on the U.S. Coast Guard’s 
Icebreaking Tug (WTGB) class cutter was studied to determine heat transfer performance 
characteristics of existing heat exchangers when used with circulating hot water vice steam. 
Characterizations such as Reynolds number vs. Colburn j factor plots, convection coefficients, overall 
coefficients, and Effectiveness-NTU relations were generated. Initial analysis with acknowledged 
conservative definitions of air side convection coefficients determined that the hydronic system 
provided on average seventy percent of the heat transfer capabilities available with the steam system. 
Improvements to the hydronic system were shown to increase heat exchanger performance 
parameters by an average of ten percent. It was notable that the added heat transfer available from 
steam is not due to a property of steam itself such as latent phase change effects, but is due solely to 
the increase in entering tube side temperature. Judging by heat transfer capabilities alone, with the 

described conservative assumptions on which these results are based, use of currently installed heat 

exchangers in a hydronic system is a viable option. 




























14. SUBJECT TERMS 15. NUMBE R OF 
Heat Transfer PAGES 156 


PRICE CODE 


17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION OF 
CLASSIFICATION OF CLASSIFICATION OF CLASSIFICATION OF ABSTRACT 
REPORT THIS PAGE ABSTRACT UL 
Unclassified Unclassified Unclassified 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 298-102 























ll 




















Approved for public release; distribution is unlimited 


ANALYSIS OF STEAM AND HYDRONIC COMPARTMENT HEATING 
SYSTEMS ABOARD U.S. COAST GUARD 140 FOOT WTGB CLASS CUTTERS 


James T. Hurley 
Lieutenant, United States Coast Guard 
B.S., Western New England College, 1984 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
June, 1996 


Author: 





Approved by: | 
Ashok Gopinath.4neais Advisor 








erry (kK. MicNelley, Chairman = \ , 
Departmertt.of Mechanical Engineering 


iii 


1V 














ABSTRACT 


The compartment heating system on the U.S. Coast Guard’s Icebreaking Tug 
(WTGB) class cutter was studied to determine heat transfer performance characteristics 
of existing heat exchangers when used with circulating hot water vice steam. 
Characterizations such as Reynolds number vs. Colburn / factor plots, convection 
coefficients, overall coefficients, and Effectiveness-NTU relations were generated. Initial 


analysis with acknowledged conservative definitions of air side convection coefficients 
determined that the hydronic system provided on average seventy percent of the heat 
transfer capabilities available with the steam system. Improvements to the hydronic 
system were shown to increase heat exchanger performance parameters by an average of 
ten percent. It was notable that the added heat transfer available from steam is not due to 
a property of steam itself such as latent phase change effects, but is due solely to the 
increase in entering tube side temperature. Judging by heat transfer capabilities alone, 
with the described conservative assumptions on which these results are based, use of 


currently installed heat exchangers in a hydronic system is a viable option. 
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I. INTRODUCTION 


A. CUTTER OPERATIONS 


The United States Coast Guard operates nine [cebreaking Tug (WTGB) class 
cutters. Stationed in the northeast United States and on the Great Lakes, these cutters are 
primarily used for domestic icebreaking, but also routinely perform other missions such 
as search and rescue, pollution response, law enforcement, and aids to navigation support. 
The cutters are 140 feet long, have a beam of 37.5 feet, displace 662 tons, ad are 
typically crewed by 17 personnel (see Appendix A for cutter illustration). The twin 
diesel-electric propulsion plant has a cruising range of 4,000 miles, maximum speed of 
14.7 knots, and provides sufficient power for the hull to break through 18 to 20 inches of 
ice. A very capable and versatile platform, these cutters and crews make vital 


contributions to the overall service the Coast Guard provides to the public. 
B. COMPARTMENT HEATING SYSTEM 


1. Purpose 

The primary purpose of the compartment heating system is to uphold the cutter’s 
overall mission readiness, keeping on board personnel physically fit and mentally alert by 
providing an atmosphere of suitable air for breathing under conditions that will enable the 
body to maintain a proper heat balance. A secondary purpose of the system is provide 
suitable temperature and humidity conditions for the preservation of stores and 


equipment. [Ref. 1] 





Zi Description of Present Configuration 

Compartment heating is accomplished either by duct heaters in the supply air 
ventilation system or by unit heaters mounted in the compartments themselves. These 
duct and unit heaters are either steam or electric. The scope of this thesis includes study 
of the steam duct and unit heaters only. 

Heating of ventilation supply air is accomplished in two stages; first by a 
preheater and then by a reheater. Preheaters are usually located near supply air ventilation 
inlets, and heat incoming air sufficiently to prevent condensation in ventilation ducts. 
Reheaters are located in the compartments being heated, and further heat the air to a set 
room temperature. Unit heaters are installed in compartments which require a spot source 
of heat. [Ref. 1] 

3; Deficiencies of Present Configuration 

The present steam system in recent years has become very maintenance intensive, 
particularly with the auxiliary boilers. The current boilers are over 15 years old, are 
“water tube” type, and are difficult and expensive to procure spare parts for. With many 
cutters operating in demanding harsh winter conditions, reduction of heating equipment 
down time is a priority. The steam heating system equipment is also non-standard 
compared to other ship classes in the fleet, leading to support and maintenance 


difficulties. 
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4. Considerations for a Replacement System 


The search for a suitable replacement for the boilers has led to consideration of a 
hydronic (circulating water in a closed piping system) compartment heating system to 
replace the current steam system. Major factors in choosing a replacement system are 
ease of installation and affordability. A way to ease the installation and decrease the cost 
of a new system is to leave in place and utilize as much of the equipment from the old 
system as possible. There is therefore considerable motivation to determine the heat 
transfer characteristics of the presently installed heat exchangers when used with 


circulating hot water. 
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Il. HEAT EXCHANGER DESCRIPTIONS 


A. SYSTEM OVERVIEW 


The compartment heating system consists of a steam generation plant and various 
heat exchangers. Included in the steam generation plant are two auxiliary boilers, 
feedwater, chemical, and condensate tanks, and associated pumps, controls, piping, and 
valves. The boilers are rated at 620 pounds of saturated steam per hour at a working 
pressure of 35 psi with feedwater entering at 210 degrees F. A functional diagram of the 
steam heating system is shown in Figure 1. 


B. HEAT EXCHANGERS 


l. General 
The heat exchangers installed on board the WTGB cutter class fall into two 
categories: unit heaters and duct heaters. Unit heaters consist of a coil, fan, motor, and 
casing assembly. They are installed in machinery and work spaces throughout the cutter. 
Duct heaters, as implied, are coils installed within the supply ventilation ducting, with 
separate ventilation fans installed upstream of coils. Duct heaters provide heat to the 
cutter’s living spaces. 
2. Unit Heaters 
a. Dimensions 
Unit heaters are manufactured by the New York Blower Company of 


Willowbrook, Illinois. Two models are employed: Model B-25 (quantity two) and Model 
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Figure 1. Steam Heating System Functional Diagram [Ref. 1] 





B-70 (quantity five). Dimensions and material properties of unit heater coils and 
assemblies are shown in Appendix B. [Refs. 2-3] 

b. Heating Capacities 

Unit heaters are sized to operate with steam at 240 degrees F and 25 psig. 


For each unit heater, the location, air flow rate, and heat transfer rate stated on the ship’s 










prints are shown in Table 1. [Refs. 4-7] 
Heat Transfer 
Stated on Prints 


(CFM) (BTU/hr) 


Table 1. Unit Heater Operating Parameters 









3: Duct Heaters 
a. Dimensions 
Duct heaters are manufactured by two manufacturers: Colmac Coil 


Manufacturing Incorporated of Colville, Washington (quantity four) and Carrier 








Corporation of Farmington, Connecticut (quantity one). Duct heater dimensions and 
material properties are shown in Appendix B. [Refs. 8-9] 

b. Heating Capacities 

Duct heaters are also sized to operate with steam at 240 degrees F and 25 
psig. For each duct heater, the supply ventilation system, location, air flow rate, heat 
transfer rate, and air temperature rise stated on the ship’s prints are shown in Table 2. 


[Refs. 4-7] 


Heat Transfer | Air Temp. Rise 
Heater System Location Air Flow | Stated on Prints | Stated on Prints 
| (CFM) (BTU/hr) (° 


01-24-1 | S01-24-1 3,680 51,600 62 to 75 


45 10 98 


Table 2. Duct Heater Operating Parameters 





C. HEAT EXCHANGER DESIGN CONSIDERATIONS 


1. Conventional Heat Exchangers 

The design of a heat exchanger involves consideration of both the heat transfer 
rates between the fluids and the mechanical pumping power expended to overcome fluid 
friction. First examining heat transfer rates, there is a marked difference in heat 


exchanger performance depending on whether the fluids involved are gases or liquids. 





The convection heat transfer coefficient, h, for gases is generally one or two orders of 
magnitude less than that for liquids. Correspondingly a gas heat exchanger’s convective 
thermal resistance, R,, defined as the inverse of the product of h and heat exchanger 
surface area, A,, (i.e. R, = 1/hA,) is one or two orders of magnitude greater than that for a 
liquid heat exchanger. It is therefore evident that for the heat transfer rate of a gas heat 
exchanger to be equivalent to that of a liquid heat exchanger, the heat transfer surface 
area for a gas heat exchanger needs to be much larger than the heat transfer surface area 
for a liquid heat exchanger. [Refs. 10-11] 

Examining the power required to overcome fluid - heat exchanger friction, 
conventional (e.g. concentric tube or shell and tube) heat exchangers operating with high 
density fluids have lesser frictional losses compared to heat exchangers operating with 
low density fluids. The pumping power required to move high density fluid (e.g. liquid) 
over a given heat exchanger is considerably less than the pumping power required to 
move low density fluid over the same heat exchanger. There can therefore be an 
equipment operating cost benefit to using a liquid heat exchanger system rather than a gas 
heat exchanger system. 

Heat exchangers where at least one fluid is required to be a gas therefore need an 
improved design to make up for inherent drawbacks. The heat transfer obtained per unit 
of heat exchanger surface area can be increased by increasing the fluid flow velocity. 
This however is not a desirable method to increase heat transfer since the friction power 
expended to increase fluid flow velocity increases by as much as the cube of velocity. 


Minimizing friction power leads to limiting flow velocities, and this combined with the 





relatively low thermal conductivity of most gases, results in low heat transfer rate per unit 
of heat exchanger surface area. There is therefore, in a conventional heat exchanger using 
gases, poor heat transfer performance at low flow velocities and an uneconomical, 

less than acceptable trade-off for increasing heat transfer per unit surface area by 
increasing flow velocities. 

2 Compact Heat Exchangers 

The development of compact heat exchangers is in response to the need to attain 
higher heat transfer rates with minimum space and power requirements. Large but 
compact surface areas are a typical characteristic of gas heat exchangers. Heat 
exchangers used on the WTGB cutter class heating system are one style of compact heat 
exchanger that incorporates large gas-side surface areas with dense arrays of continuous 
fins. 

It is first noted that compactness itself leads to high performance. A compact 
surface has small flow passages and the heat transfer coefficient, #, varies as a negative 
power of the flow passage size. A customary expression for the size of a non-circular 
flow passage is the hydraulic diameter, D,, equaling four times the cross-sectional area 
divided by the wetted perimeter. It is also true that a smaller hydraulic diameter increases 
friction power, but the benefits that compactness has on the heat transfer coefficient 
generally outweigh the detrimental influence of small hydraulic diameter on friction 
power. [Ref. 10] 

In addition to the influence of small hydraulic diameter, increases in heat 


exchanger performance can be obtained by any modification of the surface geometry that 


10 








results in a higher heat transfer coefficient at a given flow velocity. One widely accepted 
modification is use of extended surfaces or fins so that in addition to providing increased 
heat transfer surface area, the interrupted surface prevents thickening boundary layers 
from reducing heat transfer. Finned surfaces also increase friction power and thermal 
resistance due to conduction, but a small improvement in the heat transfer coefficient can 
more than offset these negative factors. [Ref. 10] 

Other methods of obtaining increased performance by change of flow surface 
geometry include the use of curved, corrugated, or wavy passages, in which boundary 
layer separation and turbulence (promoters of heat transfer) are induced. Such surfaces 
are incorporated in the duct heaters, but not in the unit heaters being studied. 

A common descriptor of compact heat exchangers is the area density, a, which is 
the ratio of heat transfer surface area to heat exchanger volume. A conventional cutoff 
for labeling a heat exchanger as compact is an area density value greater than 700 square 
meters per cubic meter (or 213 square feet per cubic feet). This is not a staunch rule 
however, as many heat exchangers have been grouped into the compact category with 
lesser area densities [Ref. 10]. 

Compact heat exchanger designs provide the benefits of high heat transfer rate 
with minimum volume and thus are very well suited for duct heater applications. In 
shipboard systems where volume savings are invariably sought, use of compact heat 


exchangers is very common. 
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Il, HEAT EXCHANGER ANALYSIS 


A. STANDARD PRESENTATION OF PERFORMANCE DATA 


AND THE REYNOLDS ANALOGY 


In engineering practice, it is often desirable to use a common presentation of 


performance data so as to avoid confusion associated with many arbitrarily defined 





parameters. In the study of heat transfer and flow-friction, commonality in data 
presentation is found using the Reynolds Analogy. As presented by Incropera and Dewitt 
[Ref. 11], and also Kays and Crawford [Ref. 12], with certain restrictions (noted shortly), 
relations that govern velocity boundary layer behavior are the same as those that govern 
the thermal boundary layer. From this it is known that non-dimensional friction and heat 
transfer relations for a particular geometry are closely related. Specifically, the Reynolds 
number, Re, the velocity boundary layer’s friction coefficient, C,, and the heat transfer 


boundary layer’s Nusselt number, Nu, are related as follows: 


Re 
C,—+ = Nu, (1) where: Re, = 2 (2) 
2 y 
hL 
Nu, = es (3) 
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Replacing Nu by the Stanton number, St and introducing the Prandtl number, Pr : 


Nu, h 
= (4) where: St 
Re, Pr p Ve, 


tH 








(5) 
cu 

Pr = 6 

r A (6) 


The relation, assuming the pressure drop in the flow direction is zero and that Pr = 1 (true 
for most gases), now takes the form: 


<f = St (7) 

This expression, known as the Reynolds analogy, relates key parameters of the 
velocity and thermal boundary layers. The accuracy of this expression depends on the 
noted restrictions, that the pressure gradient in the flow direction is zero and the Prandtl 
number equals one. It has been shown that the analogy may be applied over a wide range 
of Prandtl numbers if a correction is added. With this correction arises the modified 


Reynolds, or Chilton-Colburn analogy, shown as follows: 


C 
= = St Pr? = j, 06<Pr<60 (8) 


where j,,18 the Colburn j factor. For laminar flow, the modified Reynolds analogy is 
again only appropriate when the pressure drop in the flow direction is zero. In turbulent 
flow, conditions are less sensitive to the effect of pressure gradients and the equation 


remains valid for small pressure drops. 
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The benefits of this analogy lie in the ability to deduce heat transfer information 
from skin friction information and vice versa. A wide variety of compact heat exchanger 
performance data was compiled by Kays and London [Ref. 10] by way of plots of 
Colburn / factor and friction factor versus Reynolds number. The scope of this thesis 
includes presentation of Colburn j factor versus Reynolds number and does not include 
study of friction factor. In using this standard presentation, future correlation of friction 


factors in this standardized form is possible. 
B. MANUFACTURER’S DATA 


Plots of Colburn j factor versus Reynolds number, or any related information; 
were sought from the each of the heat exchanger manufacturers with no success. It was 
evident that in order to see desired heat exchanger performance characteristics, more 
commonly available information such as that shown in Figure 2 would have to be recast 


to the accepted non-dimensional j,, and Re forms. 


Air Temperature Rise at 5 PSIG, 0° EDB 
(WF or WR Alum. Fins) 
Face Velocity, SFPM 








Figure 2. Typical Manufacturer’s Heat Exchanger Performance Data 
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C. ANALYSIS MODEL 


The heat exchangers studied are represented in Figure 3: 





Figure 3. Heat Exchanger Representative Model ™t 


where: T, = Air (cold fluid) inlet temperature 
T, = Air (cold fluid) outlet temperature 
T; = Steam (hot fluid) inlet temperature 
T, = Steam (hot fluid) outlet temperature 


D. COMPUTATIONS AND RESULTS 


1. Relation Definitions 
Relations used in the steam system heat exchanger computations are defined as 


follows: 


Mean air temperature: 
_ (1, +T,) 
2 


Dy 


(9) 
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Air mass flow rate: 


*A *V 


M air = Pair ir air 


Aur side heat transfer rate: 


Q air ~ Mn oir ei. *(T, -f,) 


Air side convection coefficient: 


ae Q vir 
(T's 7 P) *A, 











Mass velocity: 
G at air 
Ay 
Reynolds number: 
G*D, 
Re. 2 
m 
Prandtl number: 
c *L 
Pr = 
k 
Stanton number: 
St = h 
G *c 
P 


Colburn j factor: 


jy = St* Pr? 
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(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 





2 Sequence of Computations 

Computations for the steam system heat exchanger analyses were performed with 
assistance of Fortran computer codes. The sequence of computations in the Fortran code 
titled “hxair.f’, shown as part of Appendix C, is summarized in Figure 4. Manufacturer’s 
data for each heat exchanger similar to that shown in Figure 2 were input to the program. 
Dimensional characteristics of the heat exchangers input to the program were uhied 
from manufacturer’s drawings (see Appendix B for a tabular summary of all pertinent 
dimensions). 

3. Tabular and Graphical Results 

Program output, shown in tabular form in Appendix C, included the Colburn 7 
factor and Reynolds numbers needed for further analysis of a hydronic system in a 
subsequent chapter. Plots of j,, vs. Re for each heat exchanger studied are shown in 
Figures 5 through 11. 

4. Comparison with Established Results 

Results suidined from the steam system analysis were compared with the 
established results of Kays and London [Ref. 10]. For a similarly configured heat 
exchanger, a Colburn j factor versus Reynolds number plot is shown in Figure 12. It is 
evident that orders of magnitude of Colburn / factors at respective Reynolds numbers and 
overall trends of data for this analysis compare favorably with previously established 


results. 
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Input heat exchanger dimensional 
characteristics, air flow inlet 


temperature, air velocities, and array of 
corresponding outlet air temperatures. 


Begin DO loop to compute flow 
parameters at each flow velocity. 


Compute mean air temperature to 
determine transport properties with 
FUNCTION subroutines [Refs. 13-14]. 





Using relations defined in Section III-D- 
1, compute parameters m,,, QO Mai 


iv? 


G, Re,, 





Print results to output file. 
Continue DO loop to completion. 


Figure 4. Steam System Computations Flow Chart 
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Figure 5. Colburn j Factor vs. Reynolds Number - Unit Heater B25 
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Figure 6. Colburn Factor vs. Reynolds Number - Unit Heater B70 
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Figure 7. Colburn j Factor vs. Reynolds Number - Duct Heater 01-24-1 
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Figure 8. Colburn j Factor vs. Reynolds Number - Duct Heater 01-25-1 
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Figure 9. Colburn j Factor vs. Reynolds Number - Duct Heater 01-50-0 
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Figure 10. Colburn j Factor vs. Reynolds Number - Duct Heater 2-25-1 
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Figure 11. Colburn 7 Factor vs. Reynolds Number - Duct Heater 2-16-1 
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Finned circular tubes, surface 8.0-3/8T. 





Tube outside diameter = 0.402 in = 10.2 x 103m 

Fin pitch = 8.0 per in = 315 per m 

Flow passage hydraulic diameter, 4r, = 0.01192 ft = 3.632 x 10°3m 
Fin thickness = 0.013 in = 0.33 x 10°3m 

Free-flow area/frontal area, o = 0.534 

Heat transfer area/total volume, a = 179 ft2/ft? = 587 m2/m?3 

Fin area/total area = 0.913 

Note: Minimum free-flow area in spaces transverse to flow. 


Figure 12. Established Colburn / Factor vs. Reynolds Number Relation [After Ref. 10] 
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IV. HYDRONIC SYSTEM ANALYSIS 


A. APPROACH 

To analyze the compartment heating system as a hydronic (circulating water) 
system, the Effectiveness - NTU Method is utilized. Required for this method are 
determinations of heat exchanger overall heat transfer coefficients, which entails 
computing other parameters such as air and water convection coefficients, fin 
efficiencies, wall conduction resistances, and fouling resistances. Results from the 
previous chapter are an integral part of this analysis, and are used specifically to 


determine air side convection coefficients at specified fan operating conditions. 
B. OVERALL HEAT TRANSFER COEFFICIENT 


1. Newton’s Law of Cooling 

A standard measure of a heat exchanger’s total thermal resistance to heat transfer 
between two fluids separated by a heat exchanger structure is the overall heat transfer 
coefficient, U. In an expression analogous to Newton’s Law of Cooling, the heat transfer 
between two fluids separated by one or more thermal resistances, R,, is: 


AT 


» = — = UA AT 18 
Q TR : (18) 


The result apparent from the above expression, that 1/)'R, = UA,, is applicable to 


clean, unfinned surfaces only, but can be used as the basis for determining a heat 


exchanger’s overall heat transfer coefficient with additional factors considered such as the 
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effects of fins on the air side surfaces, the presence of fouling on both air and water sides, 
and tube wall conduction. 
Zs Overall Coefficient Elements 

a. Air Side Convection Coefficients 

The air side convection coefficients, h,,,, are determined using the 
Reynolds number vs. Colburn j Factor plots, generated as part of the steam system 
analysis in the previous chapter. The utility of presenting data in non-dimensional form is 
now very apparent, as the performance characteristics of each heat exchanger formulated 
with steam now is relevant to the hydronic analysis. 

An actual operating condition of the fan supplying air to a particular heat 
exchanger, in the form of a quantity of cubic feet per minute of air at a specified air 
temperature (taken from manufacture’s data) is used to compute the air flow’s mass flow 


rate, m_., then mass velocity, G, then Reynolds number, Re. Entering the Re vs. j,, graph 


air? 


with the computed Reynolds number, a value for j,; is obtained, from which h,,, is 
computed using the Stanton number and Colburn j factor relations, simplified into the 


following: 


JyGe 
= =? (19) 





Computations of the air side convection coefficient for each of the 
system’s heat exchangers are shown in Table 3. Note that model B-25 unit heaters (i.e. 2- 


8-0 and 2-73-1) and model B-70 unit heaters (2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2) 
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each have one type fan and thus a common volumetric flow rate of air (“CFM of Air @ 


70°F”), and were each grouped into one row in the table. 


Number @70°F | (tbm/hr) | (lbm/ft’-sec) (BTU/hr-ft’-R) 
TypeB.25" | 400 | 1.798 | 116 | 2981 | oooss | 60 
‘TypeB-70" | 1505 | 6763 | 225 | 5,795 | 0.0037 | 9.0 
ST A 


a roe eer tae 
216-1 | 350 | 1573 | 082 | 851 | 0.0062 


"Type B-25 includes unit heaters 2-8-0 and TB. l. 
™ Type B-70 includes unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 









Table 3. Air Side Convection Coefficient Computations 


b. Water Side Convection Coefficients 

The water side convection coefficients, h.,,, are computed based on 
internal flow relations presented by Incropera and DeWitt [Ref. 11]. Initial computations 
of water side convection coefficients are based on a manufacturer’s recommended 
maximum water flow rate of ten gallon per minute. These initial water side convection 
coefficients will be modified in some cases based on an optimized water mass flow rate, 
as will be shown in a later section. It is also noted, as it eas in the previous chapter, that 


the water side convection coefficient is generally one or two orders of magnitude higher 
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than the gas side convection coefficient. The water side convection thermal resistance (R, 
= 1/hA) will therefore be one or two orders of magnitude lower than that of the air side, 
and a relatively small contributor to the sum of thermal resistances used to compute the 
overall heat transfer coefficient, U, as will be shown 1n a following section. 

(1) Reynolds Number for Internal Flow. The Reynolds number for 
internal flows can be computed using the following relation, where the hydraulic 
diameter, D, = 4A /P. 


4m wir 
Re, = ——— (20) 
TD, wu 


wtr 


(2) Nusselt Number. The Nusselt Number, Nup, relations for 


internal laminar flow (Re, < 2,300) and internal turbulent flow 


(Re, > 2,300 - Gnielinski Correlation) are then utilized to solve for h,,,,: 





N sree D, (21) 
u — 
° k 
= 4,36 Re, < 2300, circular duct heater tubes 
= 6.49 Re, < 2300, rectangular unit heater tubes 


(778) (Re, — 1000) Pr 
= Re, > 2300, duct heater and unit heater tub. 
1 + 12.7 (778)? (Pr? - 1) 


where: f = (0.79InRe, - 1.64)? 
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(3) Sequence of Computations. Computations of the water side 
convection coefficients for each of the system’s heat exchangers were completed with 
assistance of Fortran computer codes. The sequence of computations in the Fortran code 
titled “hxwitr.f’, shown as part of Appendix C, is summarized in Figure 13. Dimensional 
characteristics of the heat exchangers input to the program were determined from 
manufacturer’s drawings (see Appendix B for a tabular summary of all pertinent 
dimensions). 

(4) Tabular Results. Excerpts from program output is shown in 
Tables 4 through 6. Note that the water tubing dimensions for type B-25 and B-70 unit 
heaters are identical, and thus these results are grouped into Table 4. Similarly, tubing 
dimensions for all duct heaters with the exception of heater 01-24-1 are identical and 
these results are grouped into Table 5. Results for heater 01-24-1, with its own unique 
tubing dimensions are shown in Table 6. Complete program output is shown in 
Appendix C. 

c. Extended Surfaces (Fins) 
As discussed in the previous chapter, extended surfaces or fins are often 


added to heat exchanger surfaces to decrease thermal resistance to convection heat trans- 





fer. A representative fin array for the heat exchangers being analyzed is shown in Figure 


14. Expressions for fin efficiencies are as presented by Incropera and DeWitt [Ref. 11]. 
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Input tube dimensions considering tube 
may not be circular. 
Compute hydraulic diameter of tube. 


Begin DO loop to compute flow 
parameters at mass flow rates ranging 
from 0 to 4,876 Ibm/hr (10 gpm) in 
increments of 48.76 lbm/hr. 


Using relations defined in Section I[V-B- 
2-b, compute flow parameters Re,, f 
Nu,, and h 





Print results to output file. 


Continue DO loop to completion. 


Figure 13. Water Side Convection Coefficient Computations Flow Chart 




















Flow Rates 
Mass Volume 


Nup 
(lbm/hr gal/min 








828.9 1.700 25209.4 113.04 929.6 
4876.0 10.00 148290.7 492.92 4053.2 


Table 4. Unit Heater Water Side Convection Coefficient Computations 


Flow Rates 
Mass h, 
om/h gal/min 


128.6 


202.4 
908.0 


955.3 
4165.2 





Table 5. Duct Heater Water Side Convection Coefficient Computations (with exception 
of heater 01-24-1) 
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Flow Rates 
Mass Volume 


(gal/min) (BTU/hr-ft’-R) 


(Ibm/hr) 


Table 6. Duct Heater 01-24-1 Water Side Convection Coefficient Computations 















(1) Single Fin Efficiency. The single fin efficiency, n,, is given by the expression: 





tanh mL on 
— (22) 
ML ig, 
where: 
hi? 
m= (23) 
Kin A ep, 
tin 
Lin = Ln + “ (24) 
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Figure 14. Representative Fin Array 


(2) Fin Overall Surface Efficiency. Characterizing the 
performance of an array of fins is done with the expression fin overall surface efficiency, 


n,, where: 


q q 
nN, = ns (25) 
a hA, (T ) 


aie base oo 


Breaking down the total heat transfer into contributions from the 


finned and unfinned (base) surfaces and noting that A, = A, +A, .,> this expression can 
be further simplified yielding: 
A, 
ee ame ~N,) (26) 


S cir 


(3) Sequence of Computations. Computations of single and 
overall fin efficiencies for each of the system’s heat exchangers were completed with 


assistance of Fortran computer codes. The sequence of computations in the Fortran code 
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titled “hxU.f’, shown as part of Appendix C, 1s summarized in Figure 15. Dimensional 
characteristics of the heat exchangers input to the program were determined from 
manufacturer’s drawings (see Appendix B for a tabular summary of all pertinent 
dimensions). 

(4) Tabular Results. A summary of program output for each of the 
system’s heat exchangers is shown in Table 7. Efficiencies were computed using 
operating conditions and air side convection coefficients of Table 3. Complete program 


output is shown in Appendix C. 


[ewer 


B-70° Unit Heaters 0.84 





Table 7. Heater Single Fin and Overall Fin Efficiencies 


d. Fouling Resistance 
During normal heat exchanger operation, surfaces are often subject to 


fouling by fluid impurities, rust or scale formation, and other reactions between the fluid 
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and the wall material. Depositions such as these can greatly increase the resistance to 
heat transfer between the fluids. This effect can be treated by introducing an additional 
thermal resistance, termed the fouling factor, R,, the value of which depends on the 
operating temperature, fluid velocity, and length of service of the heat exchanger. 
Fouling factors recommended by the Tubular Exchanger Manufacturer’s Association are 
shown below in Table 8. The shaded fouling factors are those used for the heating system 
analyzed in this thesis. 

é. Tube Wall Conduction 


Derivations of the tube wall conduction resistance, R,,,,, by Incropera and 


wall 


Fouling Factor (hr-ft?-°F/BTU) 
Below 125 °F Above 125 °F 








Table 8. Typical Fouling Factors 
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DeWitt [Ref. 11] produce the following expressions for duct heater tubes (cylindrical) 
and unit heater tubes (modeled as plane wall). 

a 

In| — 

d. 


Re Se (27) 
wall (2TLk) 


tube 


(cylindrical tube wall conduction resistance) 


f 
R= ( ) (28) 
ds KA , tube 


(plane wall tube conduction resistance) 
3. Summation of Elements 
The summation of thermal resistances including effects of air side convection 
with fins, water side convection without fins, fouling on both air and water sides, and 


tube wall conduction yields the following: 


ee ee eee 
UA U air A Ss U er A RY 
(29) 
Pes tas 


= ———_————— + ——_————_ + + 


(nN hA A Swat * 7 hA 
(N, ate (, aie * Ee ( ey 





The calculation of the UA, product does not require designation of the air side or 


water side, i.e. UA, = U,, A, = U,, A, . However, the calculation of an overall 


coefficient, U, does depend on whether it is based on the air side or the water side surface 


area since’ UL FU if A, FA As mentioned previously, the convection 


H wir 
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coefficient for the air side is generally much smaller than that for the water side, and thus 
dominates determination of the overall coefficient. 
4, Overall Coefficient Computations and Results 

a. Sequence of Computations 

Computations of the overall coefficients for each of the system’s heat 
exchangers were completed with assistance of Fortran computer codes. The onrren of 
computations in the Fortran code titled “hxU.f’, shown as part of Appendix C, is summa- 
rized in Figure 15. Program inputs are the previously computed values of h.,,, h,,,, tabu- 


air> ~~wtr? 


lated values ofR, andR, , and dimensional and material characteristics (see Appendix 


B for a tabular summary of all pertinent dimensions). 
b. Tabular Results 
A summary of program output for each of the system’s heat exchangers is 
shown in Table 9. Complete program output is shown in Appendix C. 
C. THE EFFECTIVENESS - NTU METHOD 
1. Applicability 
The Effectiveness - NTU analysis method is used to take the overall heat transfer 
coefficient parameter and determine heat exchanger heat transfer performance 
characteristics. This method is chosen over the /og-mean temperature difference method 
based on a detailed comparison between the two methods presented by Kays and London 


[Ref. 10]. In the comparison, steps for each method are evaluated for a variety of 


4] 


Input water and air side convection 
coefficients and heat exchanger 
dimensional characteristics. Classify 
computations as either for a unit heater 
or duct heater. 


Begin DO loop to compute parameters 
for each air side convection coefficient 
input. 


Using relations defined in Section IV-B- 
2-c, compute single fin efficiency and 
overall fin efficiency. 


Using relations defined in Sections IV- 
B-3, compute overall heat transfer 
coefficient. 





Figure 15. Fin Efficiency and Overall Coefficient Computations Flow Chart 
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(BTU/hr-f?-R) 
: 4 
: 6 
3 
5 







-2 5. 
- 8. 
4. 


2-16-1 Duct Heater a ae 


" Type B-25 include unit heaters 2-8-0 and 2-73-1. 7 
™ Type B-70 include unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 


Table 9. Heat Exchanger Overall Coefficient Results 


analysis starting points. For this analysis, where only the inlet hot and cold fluid 
temperatures are known, the Effectiveness - NTU method is more straightforward, 
eliminating a tedious, iterative procedure where outlet temperatures would be estimated 
and then adjusted until the heat transfer rate corresponds to the inlet/outlet temperature 
difference. In the Effectiveness - NTU (abbreviation for “number of transfer units”’) 
approach, from knowledge of the heat exchanger type, size, and fluid flow rate, the NTU, 


maximum heat transfer rate, gQ __, and effectiveness, € (both defined shortly) are used to 
max 


determine the actual heat transfer rate. 
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2. Maximum Possible Heat Transfer Rate 

The maximum possible heat transfer rate for a given heat exchanger is achieved 
when, by definition, the entering cold air temperature is raised to equal its highest 
possible value, i.e., that of the entering hot water temperature. Referring to Figure 3, 


entering ait would rise in temperature by the quantity (T, - 7,). Using a relationship for 
the heat transfer rate similar to that presented earlier: 

O = mc,(T, - T,) (30) 
where the term mc, is defined as the heat capacity rate, C, and will be either based on 
the mass flow rate, m , and specific heat, Ca. for air or for water. As developed by 


Incropera and DeWitt [Ref. 11], the maximum heat transfer rate, Q __, is then: 


QO nnax = C min (13 7 T,) (31) 


whereC_,, is equal to eitherC,,, or C,,,, whichever is smaller. 


wtr? 


3. Effectiveness 
The effectiveness, €, is defined as the ratio of the actual heat transfer rate for a 


heat exchanger to the maximum possible heat transfer rate. 


Ee = (32) 


Q 
O sax 


The effectiveness, which is dimensionless, must be in the range 0 < € < 1. Itis 


























useful because, if . T,,and T, are known, the actual heat transfer rate may be 
determined from the expression: 
Q=e€Q . =€C_.(T, ~ T,) (33) 
4. Number of Transfer Units 


The number of transfer units, NTU, 1s a dimensionless parameter that is widely 


used for heat exchanger analysis and is defined as: 





UA, 
NTU = (34) 
C 


3: Effectiveness - NTU Relations 
The Effectiveness - NTU relationship for the type heat exchangers being analyzed, 


1.e. a single pass cross flow heat exchanger with both fluids unmixed is: 





a l 0.22 7 0.78] _ 
€ = 1-exp wre {exp |-C, (wru)>*| - 1} (35) 
where: 
C= Ci 
Sa (36) 
Cac 


The NTU, e, C, relation may be represented graphically as shown in Figure 16. 
6. Computations and Results 
a. Sequence of Computations 


Computations involving the Effectiveness - NTU analysis of the system’s 
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heat exchangers were completed utilizing a spreadsheet format, shown in Appendix D. 
Inputs are heat exchanger dimensions, water and air inlet temperatures 7, and T; 


respectively, air and water volumetric flow rates, and previously determined values of j,, 


‘ 
{ 
{ 
} 
ee eee 
i 
t 
{ 


| 


or 





Figure 16. NTU, é€, C, Relation of a Single-Pass Cross Flow Heat Exchanger with Both 
Fluids Unmixed 


Agin Ayyp and U,,. Using relations defined in Sections IV-C-2 through IV-C-5, NTU- 
effectiveness parameters were computed and finally the heat transfer rate and outlet air 
and water temperatures T, and T, were determined. 

b. Tabular Results 


Effectiveness - NTU analysis results for the hydronic system are shown in 


Table 10. 
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Water Flow Parameters 






Air Flow Parameters 


Heater 





) 2s | 400 | 1798 | 1.16 | 2981 | 0.0048 6.0 | 431 4 
ce 1505 | 6763 | 2.25 | 5795 | 0.0037 EI 


-50-0 | 750 
2-25-1 | 1050 | 4719 3654 | 0.0054 1132 | 4876 | 4053 | 4886 pa 169873 22526 60 185 
2-16-1 1573 0.0062 377 | 4876 | 4053 | 4886 54737 | 0.33 | 0.08 | 0.28 | 15152 187 


Table 10: Effectiveness - NTU Analysis Results - Hydronic System 
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V. SYSTEM COMPARISONS AND DISCUSSION 


A. OUTLINE 


The hydronic system heat exchanger performance data generated in the previous 


chapter will be compared with steam system performance parameters. Specifically, the 


heat rate outputs (Q ) and the final air temperatures (T7,) of the hydronic system will be 


compared with those of the steam system. In cases where the hydronic system 
performance matches or exceeds steam system performance parameters stated on the 
ship’s prints, water flow rates will be optimized. In instances where the hydronic system 
performance falls short of that of the steam system, further analysis will be conducted to 


attempt to increase hydronic system heat transfer. 


B. FURTHER EXAMINATION OF STEAM SYSTEM 
CAPABILITIES 


1. Modifications to Effectiveness-NTU Method 

Steam system heat exchanger heat transfer performance data is essential in order 
to assess the capabilities of the hydronic system. To compute steam system performance 
parameters, the Effectiveness-NTU analysis method is again used, modified for steam as 


follows: 





1. The convection coefficient for steam is taken to be very large (on the order of 


to the overall heat transfer coefficient. 


10*), thereby decreasing the air side thermal resistance, contributing negligibly 
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2. The specific heat for a substance undergoing a phase change, such as steam 


condensing, is infinite. The heat capacity rate for steam, C,,.,,,,. 18 therefore 


steam? 


also infinite, yielding a heat capacity ratio, C,, effectively equal to zero. 


3. In the case of C, =9 , the effectiveness, © = 1 - exp(-NTU )_, 


4. The inlet air and steam temperatures are as stated previously in Chapter II. 


5. The steam mass flow rate, M eam = O Kh»)? where h,, is the enthalpy 


difference due to condensation from saturated steam to saturated water. 
2. Sequence of Computations 
Computations involving the Effectiveness-NTU analysis of the system’s heat 
exchangers were completed utilizing a spreadsheet format shown in Appendix D. Inputs 
are heat exchanger dimensions, steam and air inlet temperatures, T, and 7’; respectively, 


as well as air volumetric flow rates, and previously determined values of j,,, h,;,. and U,,,. 


Using relations defined in Chapter III, Effectiveness-NTU parameters were computed 
and finally the heat transfer rate and outlet air temperature 7, for each heat exchanger 
were determined. 

3. Tabular Results 


Effectiveness-NTU analysis results for the steam system are shown in 


Table 11. 
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Heater 











Air Flow Parameters Steam Flow Parameters 
” Or. 
Unit NTU C, T,| 
BTU BTU 
ee aw) | | en | EP 
a i ft? q| 








B-25 
B-70 
01-24-1 
01-25-1 
01-50-0 
2-25-1 


2-16-1 


rae Tome [oxe[om [ow] oo Tor [vs Too TT so [re [one [om [ee [roma [oT 
oe [os [5 [ows [v0 [ins [as [mo [= [or [anne [es foo fee [ous [os [om | 
woo fea ise fos [oo | = [os [wos [oat [om fon [we [or [a 
| 750 3371 1.16 0.0056 809 59.0 P50 | ow 

Pe 

a 


1050 | 4719 | 1.73 | 3654 | 0.0054 | 10.1 | 1132 | 33.0 | >>0 
350 _ 1573 | 0.82 0.0062 377 | 22.5 | >>0 


Table 11: Effectiveness-NTU Analysis Results - Steam System 





C. COMPARISONS OF HYDRONIC AND STEAM SYSTEMS 


1. Tabular Results 





Comparisons of the hydronic and steam systems’ heat transfer performances are 
shown in Tables 12 and 13. 
Q hydronic Q hydronic 


O steam 


11,713 16,722 
28,214 40,445 
O siag 


Table 12: Unit Heater Hydronic vs. Steam Performance Data 
O AT 
Heater | (BTus/hr) (BTU) | (°F) | 
01-24-1 74,778 
2 


orsno fas | ss) sass | | an | am 
57 | | 























Q hydronic 


GIP pydronic 


AT 


air 


















O ie 







21428 | s7__| 070 | 070 __ 


Table 13: Duct Heater Hydronic vs. Steam Performance Data 
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2. Comparison Trends 

Comparisons between hydronic and steam system results bear some consistencies. 
Notably, the ratio of the hydronic to steam heat performance values average 
approximately 70 percent. This is a performance difference which warrants examination 
to determine what the possibilities are for improved performances and what are the causes 


for the shortfalls. 
D. IMPROVING PERFORMANCE OF HYDRONIC SYSTEM 


1. Altering Air Mass Flow Rate 

a. Analysis Method 

The Effectiveness-NTU analysis makes very apparent that the air-water 
heat exchangers studied are “air side limited”. This can be explained qualitatively by 
examining, for given heat exchanger dimensions and flow conditions, the relative 
magnitudes of the contributors to the basic relations used in the analysis. The inverse 
relationships of thermal resistances in the computation of overall coefficient, U, heavily 
skews the resulting value of U to the lower convection coefficient in the equation , in this 
case h,,,. The number of transfer units, NTU, is directly proportional to U, thus a lower 
value of U results in a lower NTU. Examining the graphical representation of the N7U-e- 
C, relation (see Figure 16), lower values of NTU yield lower values of effectiveness €. 


The heat transfer rate, , is then the product of g _ and €. 


A means to increase the air side convection coefficient, h,,, and thereby 


air? 


increasing the overall coefficient and the heat transfer rate, is to increase the rate of air 
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flow passing over the heat exchanger. A change in the air side convection coefficient 
with increasing air flow rate can be predicted by examining the plots of Colburn j factor 
vs. Reynolds number and the relation from which h,,,, is derived, i.e. 


hay = Gy Ge, (Pr 2/3) ‘The negative slope of the plots indicate that j,, decreases with 


increasing Reynolds number, but it is also true that the mass velocity, G, increases with 
increasing Reynolds number. The greatest increases in h,,,. will be realized when Re and 
therefore G increase with minimal decrease inj, This leads to the qualitative deduction 
that increases in h,,, are more profound if the plot of Colburn / factor vs. Reynolds 
number is nearer to horizontal. 

b. Analysis Results 

The net effects on h,,, of the counteracting trends of j,, and G, and the 
effects of the plot’s slope can best be seen with sample computations. Results of 
Effectiveness-NTU analysis computations at air flow rates of up to twice the original flow 
rates for duct heater B-25 are shown in Table 14. Computations for the remaining heat 
exchangers are shown in tables in Appendix D. Plots of heat transfer rate vs. air flow rate 
for each heat exchanger are shown in Figures 17 through 23. 

Cc. Air Temperature Rise Considerations 

The increased heat transfer rate at increased air flow rate produces a 
decrease in the rise in air temperature. This can be explained by examining the 
relationship between the heat transfer rate and the change in air temperature presented in 


a previous chapter, 1.¢ T,=T, + (O/c_,)- The heat transfer rate , 9, increases but the 
air 
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Water Flow Parameters 
C 





Air Flow Parameters 






in h 






my 





(°F) ss 





hr 
71 
447] 0.0042 647 4876 | 4053 | 4886 84128 0.20 | 0.13 | 0.18 | 14807 


230 0.0035 far [aes [are [ass [oese | oe [zie [ors [oie [ons Lear [ve Lier 


Table 14: Effectiveness-NTU Analysis Results - B-25 Unit Heater with Increasing Air Flow Rate 





ag 
ON 
\O 
ON 
NT |= 
olin ita 
Gs | 7] wr 









Heat Transfer Rate (BTU/hr) 


2.9 


N 


1.5 
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Figure 17. Heat Transfer Rate vs. Air Flow Rate - Unit Heater B25 
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Figure 18. Heat Transfer Rate vs. Air Flow Rate - Unit Heater B70 
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Figure 19. Heat Transfer Rate vs. Air Flow Rate - Duct Heater 01-24-1 
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Figure 20. Heat Transfer Rate vs. Air Flow Rate - Duct Heater 01-25-1 
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Figure 21. Heat Transfer Rate vs. Air Flow Rate - Duct Heater 01-50-0 
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Figure 22. Heat Transfer Rate vs. Air Flow Rate - Duct Heater 2-25-1 
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Figure 23. Heat Transfer Rate vs. Air Flow Rate - Duct Heater 2-16-1 


62 





1000 























air’s heat capacity, C 


air? 


also increases with larger air flow rates, resulting in a net 


reduction in air temperature rise. 
The reason for providing heat to a particular volume of air needs to be 


considered to determine whether a decrease inAT,,, is acceptable at a higher ¢. If the 


purpose of providing heat to a volume of air is solely to increase the air temperature of 


potentially very cold air, as is the case with duct preheaters, then the decrease n AT, at a 
higher g would not be acceptable. If the purpose of providing heat to a space is to 


provide incident heat to surfaces in the space, with the air temperature already at 
reasonably comfortable level, (e.g. 60°F), as could be the case with unit heaters, then the 


decrease inAT,,, at a higher g would be acceptable. In the latter case, the air provided to 
the space at an increased ¢-has increased enthalpy, thereby able to increase the enthalpy 


of a volume of ‘unheated’ air to a proportionally higher value. 

d. Practicality 

Increasing the flow rate of air passing over the duct heaters, if involving 
the replacement of existing fans with larger fans, is not seemingly a practical endeavor. 
The expense, effort, and potential ramifications of duct fan replacement make this course 
of action, at an initial evaluation, not viable. 

Increasing the flow rate of air passing over the unit heaters, due to much 
better accessibility and smaller scale of the project, is seemingly a practical endeavor. An 


initial review of fans available from the unit heater manufacturer, New York Blower 
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Company, indicates that fans with larger flow rates having similar dimensions as those 
presently installed are available. An excerpt from a New York Blower technical 
publication showing various size/flow rate fans is included in Appendix E. 
2 Increasing Entering Water Temperature 

a. Analysis Method 

The hydronic system analysis performed thus far was done assuming an 
entering water temperature, 7;, of 190°F, an industry accepted hydronic system operating 
temperature. It should be noted that a non-pressurized water system will never match the 
heat transfer performance of a steam system with all other system parameters equal. This 
is due to the difference in the water’s and steam’s entering temperature, which as was 
shown previously, directly affects the Effectiveness-NTU computation of the maximum 


heat transfer rate, Q which in turn affects the heat transfer rate, ¢, (recall that 
max 


O wax ~ C min 73 . r,) and Q ~ EQ ax? 


b. Analysis Results 
The effects of increasing the inlet water temperature from 190°F to 
205°F are shown in Tables 15 and 16. 
C. Practicality 
The means to provide an increased water temperature are available from 
many commercial sources. An engineering judgement needs to be made as to whether the 


energy put into additionally heating the water is worth the energy received in the form of 
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o> =190°F Oo, =205°F o>. =205°F 


(BTU/hr) (BTU/hr) 0, sine 


11,713 13,064 


28,214 31,469 
" Type B-25 include unit heaters 2-8-0 and 2-73-1. 


™ Type B-70 include unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 
















Table 15: Unit Heater Hydronic Performance Data with Increased Entering Water Temp. 
O T,=190°F 


Heat 


: 22 
44 





QO, =205°F 


(BTU/hr) 


ITT =190°F 


(°F) 

















aur | 
16,720 | 44 


Table 16: Duct Heater Hydronic Performance Data with Increased Entering Water Temp. 





heated air. A control system that would allow water temperatures to be easily adjusted to 
higher temperatures in colder weather is desirable. To minimally impact a hot water 


heating system on the whole, local electrical booster heaters could be installed 
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near heat exchanger water inlets. This could be a desirable means to provide hotter water 


to the duct preheaters, which directly heat weather-intake air. 
E. KEY FACTORS IN HEAT EXCHANGER PERFORMANCE 


iP Dependence on Air Side Convection Coefficient 
Results of Effectiveness-NTU analysis for both the hydronic and steam system 


heat exchangers are dependent on the air side convection coefficient, h,,,.. Retracing the 


path in obtaining heat transfer rates and air temperature rises from the Effectiveness-NTU 
analysis, it was discerned that the most influential characteristic of a given heat exchanger 
was the air side convection coefficient. Due to the its relatively small magnitude, the 
value of h,,, dominated computation of the overall coefficient, U,,,. From U,,,, was 


determined the number of transfer units, NTU. Larger values of NTU yielded larger 


efficiencies, €, which when multiplied by the maximum heat transfer rate, O ox , ylelded 


increased heat transfer rates, Q. 


Gas-liquid heat exchanger performance is thus regarded as air side limited, with 
the air side convection coefficient the key parameter. The manner in which h,,,, is 
computed is worth revisiting to determine influencing factors. Recalling from Chapter 


Hl: 


: Q 
= 37 
ww (T,-T,) *4,. Se 


iv 
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In this expression, the use of T, (inlet air temperature, the lowest temperature in 
the system) results in the largest difference when subtracted from 7,, thus giving the most 


conservative (lowest) values of h,,.. Carrying this conservative value of h,,, through the 


air 
analysis as described above results in conservative values for heat transfer performance. 


2 Manufacturer’s Stated Heat Transfer Performance 


Using a manufacturer’s air temperature rise table, such as the on shown in Figure 


2, the steam system heat transfer rate can be solved for directly (i.e. Q =m ce {T,-T,] ). 


Proceeding to solve for €, NTU, U. 


air? 


and h,,, (neglecting all thermal resistances but air 


side convective resistance for simplification, i.e. 1/UA = 1/(),44),,,) will give an 


indication of what 7, equivalent a manufacturer’s h,,,, and resulting heat transfer rates/air 
temperature rises are based on. 

Looking at manufacturer’s information presented in Figure 2, and considering the 
case of a single row, 12 fin per inch heat exchanger, with an air flow velocity of 500 
FPM, the temperature rise from 0°F is 77.2°F. Solving for h,,, and then the 7, equivalent 
as outlined above yields a value of T, equal to 58°F (see Appendix F for computation 
details). This is in contrast to the value of 0°F used for computations of h,,, in the 
analysis done in Chapter III (see output for Fortran computer code “hxair” I Appendix B). 
This demonstrates that values of ,,, that manufacturer’s heat transfer rates are based on 
are significantly higher than values that were used in this thesis, leading to differences in 


heat exchanger performance values. 
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The difference in the analysis starting point of 7, outlined above indicates that 
magnitudes of heat exchanger performance data found in this thesis are conservative by 
approximately 15-20 percent when compared to manufacture’s data. The percent 
difference between the hydronic and steam system heat transfer performance will 
however still be consistent with the results of this thesis - compared data will have larger 
magnitudes but with approximately the same percent differences. The only true 
verification of any prediction of the complex heat transfer involved in a compact gas- 
liquid heat exchanger is through controlled measurements under known conditions. By 
experimental verification, predictions to cover similar geometries are then able to be 
made with more certainty. 


F. OPTIMIZATION OF WATER FLOW RATES 


iF Optimization Candidates 

Heat exchangers that meet or exceed the heat transfer rates air temperature rises as 
stated in the ship’s prints (heat exchangers 2-8-0, 2-80-1, and 01-24-1) can be further 
analyzed to determine an optimal water flow rate. With the air flow rate in this view of 
system operation fixed at the manufacturer’s/print values, the water mass flow rate can be 
adjusted to lesser values to obtain an optimal value at which the minimum heat delivery 
rate and air temperature rise can be achieved. The objective in this optimization is to 
achieve the minimum water pumping power required to meet stated system heat transfer 


performance parameters. 
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Zz Effectiveness-NTU Analysis Method Readapted 


The Effectiveness-NTU analysis method can be used once again with the intention 
of solving for the heat exchanger surface area for a given heat transfer rate. The heat 


exchanger surface area can be solved for using a relation presented in Chapter III, i.e. 





(NTU)(C i.) 
: wl min 3 8) 
™ O oir 
where, for instances of (m io (m Cy air? 
C nin ~ (m a ere (39) 


Thus, for a given heat exchanger with known air flow rate, air side convection 
coefficient, estimated water side coefficient, and stated heat transfer rate, the air side 
surface area can be computed and compared with the known value. This optimization 
procedure is continued until the computed surface area is less than or equal to the known 
surface area. 

Once an optimal water mass flow rate is found, the accompanying water side 
convection coefficient is adjusted to match the value used in the computation of the 


overall coefficient. This is not too onerous of a task since as was described in the water 





side convection computation section of Chapter III (results shown in Appendix C), the 


| water side convection coefficients at low water flow rates (i.e. 0.5 gpm) are already an 





order of magnitude above the air side convection coefficients and thus have little effect 


on the overall coefficient. 
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a. Sequence of Computations 

Computations involving the optimizations of water flow rates using the 
Effectiveness-NTU analysis were completed with assistance of Fortran computer codes. 
The sequence of computations in the Fortran code titled “hx NTUOPT.f’, shown as part 
of Appendix C, is summarized in Figure 24. Program inputs are the previously 


computed estimated values of U,,,, water inlet and outlet temperatures 7, and T, 


ir? 


respectively, air volumetric flow rates, heat exchanger air side surface areas 4, , and the 
air 


stated design heat rates 9 
Pp 


rints 


b. Tabular Results 

An example of Effectiveness-NTU water mass flow rate optimization 
program output for heater 2-80-1 is shown in Table 17. Program outputs for other 
optimized heat exchangers 2-8-0 and 01-24-1 are shown in Appendix C. The shaded last 
row represents the final optimized parameters. 

A summary of optimized water mass flow rates for unit heaters 2-8-0, 2- 
80-1, and duct heater 01-24-1 are shown in Table 18. — 

Cc. Low Optimal Flow Rates 

As shown in Table 18, for heat exchangers 2-8-0 and 2-80-1, the minimum 
water mass flow rate are very low compared to the value of 10 gpm used for other heat 
exchangers in the system. It is apparent that these two heat exchangers could have much 


less heat transfer surface area on both the water and tube sides at the expense of an 
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Input water and air inlet temperatures, overall 
coefficient, air flow rate, air side surface area, and 


heat transfer rate. 


Compute air mass flow rate and heat capacity rate. 


Begin DO loop to compute Effectiveness-NTU 
parameters at water mass flow rates iterated up to 
4.876 lbm/hr (10 gpm). 





















Compute water heat capacity rate C,,,,, and compare 
with air heat capacity rate C,,. If C,;,<C,,,. proceed 
with Ci, = Cai If Cy. < Cai, proceed with C 


min 






ad 


Compute heat exchanger air side surface area and 
compare with actual air side surface area. If 
computed value is less than actual value, print results 
to output file. 














If computed air side surface area is greater than 
actual, continue DO loop until solution is found or 
end program at water mass flow rate of 4,876 lbm/hr 


(10 gpm). 


Figure 24. Effectiveness-NTU Computations Flow Chart 
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increased water mass flow rate and still provide heat transfer rates stated on the prints. 
This can be seen from the results in Table 17, where the increase in water mass flow rate 
corresponds to a decrease in air side surface area. Water flow rates as low as those 
obtained for heat exchangers 2-8-0 and 2-80-1 delivering heat transfer rates stated on the 
ship’s prints indicates that the heat exchangers are oversized for the application. This 
apparent oversizing is also evident when it is noted that the same size unit heater (B-70) 
is used to deliver heat transfer rates varying from 5,294 to 64,860 BTU/hr. Selections of 
apparently oversized heat exchangers for a particular heat load was most likely based on 


other sound reasoning, not solely on heat transfer capacity. 


C C 
min max 
(BTU/hr-R) 





om oe) 








Table 17. Results of Water Mass Flow Rate Optimization Using Effectiveness-NTU 
Analysis for Heater 2-80-1. Actual A, = 3,013 in* and Q)yingg = 7,747 BTU/hr | 
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Actual Water Flow Rates at 
A A. (computed) S A, (actual) 















Q prints S air 


ra (BTU/hr) (in’) (Ibmv/hr) (gpm) 


Table 18. Summary of Effectiveness-NTU Optimization Analysis Results 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. SYNOPSIS 


This thesis successfully analyzed heat transfer performance aspects of a set of heat 
exchangers installed on the U.S. Coast Guard’s WTGB Icebreaking Tug class cutters. 
Initial analysis with acknowledged conservative definitions of air side convection 
coefficients determined that the hydronic system provided on average seventy percent of 
the heat transfer capabilities available with the steam system. Practical improvements to 
the hydronic system were shown to increase heat exchanger performance parameters by 
an average of ten percent. It was notable that the added heat transfer available from steam 
is not due to a property of steam itself such as latent phase change effects, but is due 
solely to the increase in entering tube side temperature. Judging by heat transfer 
capabilities alone, with the described conservative assumptions on which these results are 
based, use of currently installed heat exchangers in a hydronic system is a viable option. 


B. PROPOSED HYDRONIC SYSTEM 


A hydronic system circulating water at 10 gallons per minute and 190 degrees F, 
will provide at minimum the heat transfer and air temperature rises shown in Table 10. 
These performance parameters were arrived at considering the entire thermal circuit 
involved and with conservative assumptions (when compared to manufacturer’s methods) 
for air side convective resistance determinations. Hydronic system performance 


parameters can be improved through the practical methods outlined. Considering heat 
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transfer capabilities, use of the currently installed heat exchangers as part of a hydronic 
system is plausible. 

Use of a pressurized water system, where system pressure would allow the 
circulating water temperature to equal or exceed that of steam (240°F) is another 


approach to the hydronic system design that could be explored. 
C. FURTHER STUDY 


There are numerous issues to examine further in order to take into account all 
aspects of using presently installed heat exchangers as part of a hydronic system. Issues 
remaining are related to the hardware involved as well as further analysis/design work. 

One hardware issue remaining is the internal heat exchanger component’s 
compatibilities with circulating hot water as opposed to steam. In steam heat exchangers, 
there is typically a strip of shaped metal which distributes the flow, insuring that the fluid 
passes nearly equally through all tubes. The erosion effects of water on this and other 
internal elements designed for use with steam should be examined. 

A hardware issue on the external side of the heat exchanger is the relative location 
of water vs. steam inlet and outlet piping. One manufacturer’s information indicated that 
for a steam heat exchanger, inlet piping is above the outlet to allow for condensate 
drainage. Conversely for circulating water, the inlet is below the outlet to help prevent air 
in the system from inhibiting water flow. The effort involved in possibly reversing inlet 
and outlet piping when switching from a steam to hydronic system thus needs to be 


examined. 
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Further design and analysis work includes among other matters, the hydronic 
system design to provide water at either 10 gallons per minute or the optimized values to 
the respective heat exchangers. Pressure drops across each heat exchanger and across 
other system components need to be accurately gauged to arrive at a required pumping 
head. Air side pressure drops also need to be examined if any of the heat transfer 
enhancements involving increases in air flow rate are adopted, to insure fans remain 


properly sized. 
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APPENDIX B. HEAT EXCHANGER DIMENSIONS 
AND MATERIAL PROPERTIES 
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Type B-25 See unit heaters 2-8-0 and 2-73-1. 
™ Type B-70 includes unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 


Table B-1. Heat Exchanger Dimensions and Material Properties 


Hydraulic 
Diameter 






Area 
Density 
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| 
Tube Tube Wall eee Tube Side 
Heater Thickness onductivity | Surface Area 
Material 
f tube (in) ‘ 
/ 


$8 





" Type B-25 includes unit heaters 2-8-0 and 2-73-1. 
™ Type B-70 includes unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 


Table B-2. Heat Exchanger Dimensions and Material Properties 












Fin Thermal 
Conductivity 














Fin Thickness 
te, (in) 


- i ; ; 118 
Z 18 










2-16-1 0.438 1.299 0.0065 


" Type B-25 includes unit heaters 2-8-0 and 2-73-1. 
“" Type B-70 includes unit heaters 2-80-1, 2-40-1, 2-40-2, 2-60-1, 2-60-2. 





Table B-3. Heat Exchanger Dimensions and Material Properties 
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=¢ ) COLMAC STEAM COILS 
Z “COIL 5/8” OD TUBE 
Manufacturing Inc. 
STEAM HEATING COILS 
5/8” OD TUBE 
Dimensions For Steam Coils 
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NOTE: 
1 5/16° DIA HOLES ON 3° CTRS FROM CENTERLINE OF CASING. 


2, CDIL PITCHED IN CASING’ TOWARD RETURN END 1/4% PER FOOT OF FINNED LENGTH. 
B DIM PLUS C DIM PLUS FIN HEIGHT EQUALS H. 
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Carrier 





UNIT SIZE 398. { 04000~— 050 060 


AIR QUANTITY (Ctm 
Face Velocity (fpm) 400 













Face Velocity (fpm) 800 
(H & V units only) 900 







Face Area (sq ft) . 
U-Bend 2.15 314 36 
Steam Distributing Tube 1.73 2 


| 
t 
HUMIDIFIER CONN. 
(no. ... size) 
Atomi zing Spray 1.¥ Lat Pee 4 
I 


Steam Grid 
Supply (in. OD) 1,..1% Loy, Veco d% 


Drain (in. OD) 1 ¥ ! 1 ¥ 


FACE AND BYPASS ; : : 
DAMPER AREAS (sq ft) | 6.0 7.1 | 11.0 
MIXING BOX 
DAMPER AREAS (sq ft) | 5.9 ee 6.9 | 
FANS (no. wheels...no. outlets) ini ; _ | x 
Diam** (in,) oY | 10% 12% 
| 






Wheel Length (in.) 7% 8 9% 
No. of Blades 43 48 43 
Shaft Critical Speed (rpm) 2700 2550 1700 
Max Operating Speed (rpm) 2160 2040 1360 
Fan Sheave Bore (in.) 1%, 1%. IK, 











































































































































FAN & COIL SECTION L 19% { 40% | 53% 39% 22% | 33% | 67% 49%, 28% ° | 42 
W/ 38% 38% | 42% | 38% 38%, 38% | 42% | 46% | ay 46% | 49¥ 
H 31% 1 20% | 36% | 36% | 57% | 23% | 43 | 43 70 29%, 
MIXING BOX ION LT 21% 2 21% | 2K 24Y, 24% 1 24%" | 24H) 30%, 30% 30% | 30% 
Wl 36% 36% 36% | 36% | 36% 36% | 36% | 36% | aay, 44¥, 44¥, 44¥, 
H | 25% 25% | 25% 1 25% | 20% 29% |} 29% | 29% | 37% | 37% 37% | 37% 
4 {| 4% 4%, ‘4 % A 
FILTER SECTION Wi 36% | 36% | 36% 36% | 36% | 36% | 36% | 36% | aay | 44% | 44% | ay 
A 20% | 20% | 20% | 20% | 29% | ax | my 23% | 2% | 20% | 29% | 20% 
PREHEAT SECTION L 6 aa 7 4 7 7 = 9 9 9 = 
WI 36% | 36% 36% | = | 26 36% | 28% | a 44Y, | 42% 44Y, | 2 
oe H | 20% 20% 20% : ~ 23% 23%, 23% - 29%, 29% 2A 1 ~ 
’ “BY-PASS PLENUM: LT 6 6 6 7 7 7 & 9 OA Get fe 
FACE AND BYPASS W) 36% | 36% | 36% - 36% 36% 36% “ 44y, | ay, | aay, | - 
DAMPER. SECTION H | 28% | 28% 5 = | 28% 32% 132% | es 4t% | 41% | 41% - 
SERVICE AREA-— = A} 46% 7 3a 19% | 40% | 53% 39% | 22% ° 733% 1 67% | 49% | 28% | 4 
REQUIREMENTS: °.. B} 48% | agy | 9x 48¥, | aay 48% | 48% | 48x | 55% 1 56% | 56% + 56% 
tL — Length H — Height 
W — Width H&V — Heating and Ventilating 
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APPENDIX C. FORTRAN COMPUTER CODES 
AND OUTPUTS 
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Jim Hurley 


Naval Postgraduate School 


Spring 1996 
PROGRAM HEAT EXCHANGER MODEL AIR SIDE 


REAL CPAIR, TSTD,T1,TS,AFR,AFF,ASAIR, DH, VAIRL, VAIRU, VAIR, 
+T2 (200:1200) , TM, MDOTAIR, QDOTAIR, HAIR, G, RE, PR, ST, JH, DAIR, 
+MUAIR, KAIR, PATM, R,MU0,TO,S,C(0:6),N,VAIRINC 

DATA CPAIR, TSTD/0.24,70.0/ 


Input heat exchanger dimensional characteristics, air flow inlet 
temperature, air velocities, and array of corresponding outlet 
air temperatures. 


OPEN(55, FILE = ‘RESULTS.HXAIR’ } 
WRITE (*, 4) 
WRITE(*,*)"Air flow entering temperature (degrees F):" 
READ *,T1 
WRITE (*,*) "Heat exchanger surface temperature (degrees F):" 
READ *,TS 
WRITE (*,*) "Heat exchanger frontal area (square inches) :" 
READ *,AFR 
WRITE(*,*) "Heat exchanger free flow area (Square inches):" 
READ *,AFF 
WRITE(*,*) “Total air side surface area (square inches) :" 
READ *,ASAIR 
WRITE(*,*)*Flow passage hydraulic diameter (inches):" 
READ *,DH 
WRITE(*,*) "Coil face velocity lower value (FPM):" 
READ *,VAIRL 
WRITE(*,*}"Coil face velocity upper value (FPM):" 
READ *,VAIRU 
WRITE(*,*) "Coil face velocity increment (FPM):" 
READ *,VAIRINC 
DO 20 VAIR = VAIRL, VAIRU, VAIRINC 
WRITE (*,5)T1,VAIR 
READ *,T2(VAIR) 
CONTINUE 
WRITE (55,6)T1,TS, AFR, AFF, ASAIR, DH 


Begin DO loop to compute flow parameters at each flow velocity. 


DO 21 VAIR = VAIRL, VAIRU, VAIRINC 


Compute mean air temperature to determine transport 
properties with FUNCTION subroutines. 


TM = (T2(VAIR)+T1}/2.0 


Using relations defined in Section III-D-1, compute parameters 
MDOTAIR, QDOTAIR, HAIR, G, RE, PR, ST, and JH. 


MDOTAIR = DAIR(TSTD) * (AFR/144.0) * (VAIR*60.0) 
QDOTAIR = MDOTAIR*CPAIR* (T2 (VAIR) -T1) 

HAIR = QDOTAIR/ { (TS-T1) * (ASAIR/144.0)) 

G = (MDOTAIR*144.0)/(3600.0*AFF) 


RE = G* (DH/12.0) /MUAIR(TM) 
PR = CPAIR*MUAIR (TM) /KAIR(TM) 
ST = HAIR/ (3600.0*G*CPAIR) 
JH = ST*PR**{(2.0/3.0) 


Print results to output file. 
WRITE (55,7)}VAIR,T2 (VAIR) ,MDOTAIR, QDOTAIR, HAIR, RE, JH 


Continue DO loop to completion. 


21 CONTINUE 


4 


FORMAT (//,"****HEAT EXCHANGER MODEL ~ AIR SIDE****",//, 


AAANAN 





+“Enter the following parameters in the units",/, 
+"indicated:",/) 


5 FORMAT(/,"T2 at Tl = ",F5.1," F and VAIR = ",F5.0," FPM:") 
6 FORMAT(////, 

+T25, "HEAT EXCHANGER MODEL - AIR SIDE",/, 

4725, "SsssSssssseasssscansssscusesese=",////, 


+T5, "Heat exchanger number:",////, 


+T8,"T1 = “,F5.1," degrees F",/, 

+T8,"TS = ",F5.1,*" degrees F",/, 

+T8, "AFR = ",F5.1," square inches", /, 

+T8, "AFF = ",FS5.1," square inches",/, 

+T8, "ASATR = ",F7.1," square inches", /, 

+T8, "DH = ",F4.2," inches",///////, 

+T1, "VATR T2 MDOTAIR QODOTAIR H-AIR RE JH",/, 
+T1,"(FPM) {F} {lbm/hr) (BTU/hr) (BTU/hr sqft R)",/, 

4T1, "--------- 2 2 - - ~~ - + ee ee ee "} 


7 FORMAT(T1,F5.0,T8,F5.1,T15,F6.0,T24,F7.0,T36,F6.1,T46,F6.0, 
+T55,F6.5) 


re re ee ee wre ee oe eo ee a ne en nr ee re ee we ee ee ee ee ey ee ee ee ee ee ee ee ee ee ae es 
—— oe ee ee ee ne ee ee a ee 


FUNCTION subroutines to compute temperature dependent transport 
properties of air. 


FUNCTION DAIR(TSTD} 

REAL DAIR,TSTD, PATM,R 

DATA PATM,R/14.7,53.34/ 

DAIR = (PATM*144.0)/(R* (TSTD+459.67) ) 
RETURN 

END 


FUNCTION MUAIR(TM) 

REAL MUAIR,TM,MU0,TO,S 

DATA MUO,T0,S/1.1491E-5,491.67,198.72/ 

MUAIR = MUO*({(((TM+459.67) /TO)**1.5)* ( (TO+S) /( (TM+459.67)+S) }) 
RETURN 

END 


FUNCTION KAIR(TM) 

REAL KAIR,TM,C(0:6),N 

DATA €(0),C(1),C(2),C(3),C(4),C(5),C(6)/ 

+-2.276501E-3,1.2598485E-4, -1.4815235E-7,1.73550646E-10, 

+-1.066657E-13,2.47663035E-17,0.0/ 

KAIR = 0.0 

DO 24 N= 0.0, 6.0, 1.0 

KAIR = C(N)*(((TM+459.67)*0.5556)**N) + KAIR 

24 CONTINUE 

KAIR = (KAIR*0.57782) /3600.0 

RETURN 

END . 
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HEAT EXCHANGER MODEL ~- AIR SIDE 


Heat exchanger number: B-25 (2-73-1 and 2-80-1) 


0.0 degrees F 
227.0 degrees F 
100.0 square inches 
61.7 square inches 
3013.0 square inches 
0.38 inches 


nea we on a ow 


T2 MDOTAIR QDOTAIR H-AIR 
(F} (lbm/hr) (BTU/hr) (BTU/hr 


5 4.4 
3 5.3 
3 6.0 
9 6.6 
vi 7.2 
0 Teo 
7 738 
2 8.0 
2 8.1 








HEAT EXCHANGER MODEL ~ AIR SIDE 
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Heat exchanger number: B-70 (2-8-0, 2-40-1, 2-40-2, 2-60-11, 


T1 = 0.0 degrees F 
TS = 227.0 degrees F 
AFR = 185.0 square inches 
AFF = 120.0 square inches 
ASAIR = 5198.0 square inches 
DH = 0.38 inches 
VAIR T2 MDOTAIR QDOTAIR H~AIR RE JH 
(FPM) (F) (lbm/hr) (BTU/hr) (BTU/hr sqft R) 
400. 70.3 2310 38977 4.8 2111 .00571 
500. 66.9 2888 46365 Ln 2646 .00543 
600. 63.3 3465 52644 6.4 3185 .00514 
700. 59.9 4043 58119 7.1 3726 .00487 
800. 56.7 4620 62873 27 4269 .00461 
900. 53:0 5198 66117 8.1 4817 .00431 
1000. 49.7 5775 68889 8.4 5367 .00404 
1100. 46.2 6353 70441 8.6 5921 .00376 
1200. 42.9 6930 71356 8.7 6476 .00349 


2-~60-2) 
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HEAT EXCHANGER MODEL - AIR SIDE 


Heat exchanger number: 01-24~1 


0.0 degrees F 
227.0 degrees F 
693.0 square inches 
405.0 square inches 
17700.0 square inches 
0.25 inches 


MDOTAIR QDOTAIR H-AIR 
(lbm/hr) (BTU/hr) (BTU/hr sqft R) 


129391. 
145902. 
160440. 
173732. 
185674. 
197669. 
207170. 
217607. 
225551. 
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HEAT EXCHANGER MODEL - AIR SIDE 
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Heat exchanger number: 01-25-1 


Ti = 0.0 degrees F 
TS = 227.0 degrees F 
AFR = 288.0 square inches 
AFF = 155.0 square inches 
ASAIR = 7959.0 square inches 
DH = 0.15 inches 
VAIR T2 MDOTAIR QDOTAIR H-AIR RE JH 
(FPM) (F) {ibm/hr) (BTU/hr) (BTU/hr sqft R) 
200. 114.9 1798 49587. 4.0 498 .00785 
400. 84.7 3596. 73107. 5.8 1020 .00580 
600. 69.7 S395. 90240. 7.2 1548 .00478 
800. 60.0 7193. i03575:. 8.3 2080 .00411 
1000. 53.1 8991. 114580. 9.1 2614 .00364 
1200. 48.0 10789. 124290. 9.9 3150 .00329 
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HEAT EXCHANGER MODEL - AIR SIDE 


Heat exchanger number: 01-50-0 


0.0 degrees F 
227.0 degrees F 
216.0 square inches 
116.0 square inches 
5969.0 square inches 
0.15 inches 


Wow mW at tt tt 


MDOTAIR ODOTAIR H-AIR 
(lbm/hr) (BTU/hr) (BTU/hr 


.00579 
.00477 
00410 
.00363 
.00329 














HEAT EXCHANGER MODEL - AIR SIDE 
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Heat exchanger number: 2-25-1 


T1 = 0.0 degrees F 
TS = 227.0 degrees F 
AFR = 198.0 square inches 
AFF = 109.0 square inches 
ASAIR = 2868.0 square inches 
DH = 0.31 inches 
VAIR T2 MDOTAIR QDOTAIR H-AIR RE JH 
(FPM) (F) (lbm/hr) (BTU/hr) (BTU/hr sqft R) 
200. 68.1 1236 20205 4,5 1023 .00911 
400. ais Fae 8 2472 30323 6.7 2074 .00684 
600. 42.9 3709 38185 B.4 3233 .00574 
800. ch ew 4945 44742 9.9 4195 .00505 
1000. 34.4 6181 50587 11.2 5259 .00457 
1200. 31:5 7417 56076 12.4 6325 .00422 
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HEAT EXCHANGER MODEL - AIR SIDE 


Heat exchanger number: 2-16-i1 


0.0 degrees F 
227.0 degrees F 
144.0 square inches 
77.0 square inches 
3979.0 square inches 
0.15 inches 


| | 


MDOTAIR QDOTATR H-AIR 
(lbm/hr) (BTU/hr) (BTU/hr sqft R) 
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Jim Hurley 
Naval Postgraduate School 
Spring 1996 


PROGRAM HEAT EXCHANGER MODEL WATER SIDE 


REAL PI,D,CP,MU,K,A,B,DH, PR, MDOTWTR, VDOTWTR, RE, F, NU, HWTR 


DATA PI,D,CP,MU,K 
+/3.14159,60.79,1.002,77.0E-7,0.386/ 


Input tube dimensions considering tube 


OPEN{55, FILE = ‘'RESULTS.HXWTR’ ) 
WRITE (*, 4) 


may not be circular. 


WRITE (*,*)"‘A’ dimension of tubing (inches):" 


WRITE(*,*)"Note that for a circular tube, ‘A’ 


READ *,A 


Ht 


tube diameter.” 


WRITE(*,*}"’B’ dimension of tubing (inches):" 


WRITE(*,*) "Note that for a circular tube, ‘'B’ 


READ *,B 
Compute hydraulic diameter of tube. 


IF (A .EQ. B)DH A 


tube diameter.” 


IF (A .NE. B)DH = (4.0*A*B)/((2.0*A)+(2.0*B) ) 


Begin DO loop to compute flow parameters at mass flow rates 


ranging from 0 to 4,876 lbm/hr (10 gpm) in increments of 48.76 lbm/hr. 


PR = (CP*MU*32.174*3600.0) /K 
WRITE (55,5)A,B,DH, PR 
DO 21 MDOTWTR = 0.0, 4876.0, 48.76 


Using relations defined in Section IV-B-2-b, 


RE, F, NU, and HWTR. 


VDOTWTR = MDOTWTR*7.481/ (D*60.0) 


RE = (4.0*MDOTWTR*12.0)/{3600.0*PI*DH*MU*32.174) 


IF (RE .LE. 2300.0 .AND. A .EQ. B)NU 
IF (RE .LE. 2300.0 .AND. A .NE. B)NU 
IF (RE .GT. 2300.0)THEN 

F = (0.79*LOG(RE) - 1.64)**-2.0 

NU = ((F/8.0)* (RE-1000.0)*PR)/ 

+ (1.0 + 12.7*SORT(F/8.0)* (PR** 

ENDIF 
HWTR = (NU*K*12.0) /DH 


Print results to output file. 


4.36 
6.49 


(2.0/3.0) @1..0)} 


WRITE (55, 6)MDOTWTR, VDOTWTR, RE, NU, HWTR 


Continue DO loop to completion. 


21 CONTINUE 


4 FORMAT(//,"****HEAT EXCHANGER MODEL - WATER SIDE****",//, 


+"Enter the following parameters in the 
+"indicated:",/) 

5 FORMAT(////, 
+T25, “HEAT EXCHANGER MODEL - WATER SIDE 
+725 , “BSSSSsseSessessressseeseesarsss=5 
+T5,"Heat exchanger number:",///, 
+T8,°A ",F5.3,"“ inches",/, 
+T8,"B ",F5.3," inches",/, 
+T8, "DH " F5.3¢" dnches";/, 
+T8,"PR = °,F5.3,///, 
+T4, " Flow Rates",/, 


units", /, 


see dee 
EAD a ar ar ae 


+T4," Mass Volume Reynolds Nusselt HWTR",/, 
+T4," (lbm/hr) (gal/min) Number Number (BTU/hr sqft R)",/, 
/) 


tT 4, men mn rr enn 


ee ee 


H 
‘ 


6 FORMAT(T4,F6.1,T15,F6.3,T24,F8.1,T736,F6.2,T49,F6.1) 


compute parameters 
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Heat exchanger number: 


A 
B 
DH 
PR 


Hoi i i 


Flow 
Mass 
(lbm/hr) 
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0.310 inches 
3.080 inches 
0.563 inches 
24315 


Rates 
Volume 
(gal/min) 
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HEAT EXCHANGER MODEL - WATER SIDE 


Reynolds 
Number 
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Nusselt 
Number 


B-25 (2-73-1 and 2-80-1) 


B-70 (2-8-0, 2-40-1, 2-40-2, 2-60-1, 2-60~-2) 


HWTR 
(BTU/hr sqft R) 
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72662. 
74145. 
15628, 
THEIL; 
78594. 
80077. 
81560. 
83042. 
84525. 
86008. 
87491. 
88974. 
90457. 
91940. 
93423. 
94906. 
96389. 
97872. 
99354. 
100837. 
102320. 
103803. 
105286. 
106769. 
108252. 
109735. 
111218. 
112701. 
114184. 
115666. 
117149. 
118632. 
120115. 
121598. 
123081. 
124564. 
126047. 
127530. 
129013. 
130495. 
131978. 
133461. 
134944. 
136427. 
137910. 
139393. 
140876. 
142359. 
143842. 
a45S20% 
146807. 
148290. 
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Heat exchanger number: 


Flow 


Ma 


ss 


(ibm/hr) 
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0.555 inches 
0.555 inches 
0.555 inches 
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Rates 


Volume 


(gal/min) 
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HEAT EXCHANGER MODEL - WATER SIDE 
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Reynolds 


Number 
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Nusselt 


01-25-11, 01-50-0, 


Number 


2-25-2, 


HWTR 
(BTU/hr sqft R) 
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HEAT EXCHANGER MODEL - WATER SIDE 


Heat exchanger number: 01-24-1 


1.035 inches 
1.035 inches 
1.035 inches 
2.315 


Flow Rates 
Mass Volume Reynolds Nusselt HWTR 
(lpbm/hr) (gal/min) Number (BTU/hr sqft R) 
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Using relations defined in Section IV-B-2-c, compute single fin 
Jim Hurley efficiency and overall efficiency. 
Naval Postgraduate School 
Spring 1996 ({2.0*WFIN)+(2.0*TFIN))}/12.0 
: (WFIN*TFIN) /144.0 
PROGRAM HEAT EXCHANGER OA COEFFICIENT WATER SORT ( (HAIR(L) *P) / (KF IN*AX) ) 
(LFIN+ (TFPIN/2.0})/12.0 
(TANH (M*LC) ) / (M*LC) 


REAL HAIR(20),HWTR, KFIN, LFIN, WEIN, TFIN, ASAIR, ASWTR, ASFASAIR, = 
= 1.0-ASFPASATR* (1.0-ETAF) 


+KTUBE, TTUBE, ATUBE, LTUBE, ODTUBE, IDTUBE, P, AX,M,LC, ETAF, ETAO, 
+UAIR, RFAIR, RFWTR, PI 


INTEGER FLAG,N,J,L Using relations defined in Section IV-B-3, compute overall heat 
DATA RFAIR, RFWTR, PI/0.002,0.0005,3.14159/ transfer coefficient. 
Input water and air side convection coefficients and heat exchanger IF (FLAG .EQ. 1)THEN 
dimensional characteristics. Classify computations as either for VAIR = 1.0/ 
a unit heater or a duct heater. ({({1.0/ (ETAO*HAIR(L) * (ASAIR/144.0)))+ 
{(RFAIR/ (ETAO* (ASAIR/144.0)))+ 
OPEN(55, FILE = ‘RESULTS.HXU’ ) (1.0/ (HWTR* (ASWTR/144.0)))+(RFWTR/ (ASWTR/144.0))+ 
WRITE (*,4) ((TTUBE/12.0) / (KTUBE* (ATUBE/144.0)))) 
WRITE (*,*) "Convection coefficient - water side (BTU/hr sqft R):" . * (ASATR/144.0) } 
READ *,HWTR ENDIF 
WRITE (*,*) "Quantity of convection coefficients - air side that" IF (FLAG .EQ. 2)THEN 
WRITE (*,*) "will be entered." . UAIR = 1.0/ 
READ *,N (((1.0/ (ETAO*HAIR(L) * (ASAIR/144.0)))+ 
DO 21 J = 1, N, 1 (RFAIR/ (ETAO* (ASAIR/144.0)))}4+ 
WRITE (*,5)7 (1.0/ (HWTR* (ASWTR/144.0)))+(RFWTR/ (ASWTR/144.0) }+ 
READ *,HAIR(d) (LOG (ODTUBE/ IDTUBE) / (2.0*PI* (LTUBE/12.0) *KTUBE) ) ) 
CONTINUE * (ASAIR/144.0)) 
WRITE (*,*) "Conduction coefficient - fin (BTU/hr ft R):” ENDIF 
READ *,KFIN 
WRITE(*,*)"Fin length (inches) :" Print results to output file. 
READ *,LFIN 
WRITE (*,*) "Fin width (inches):" WRITE (55,8)HAIR(L) , ETAF, ETAO, UAIR 
READ *,WFIN 
WRITE(*,*)} "Fin thickness (inches):" Continue DO loop to completion. 
READ *,TFIN 
WRITE(*,*) "Total air side surface area {sq. in.):" CONTINUE 
READ *,ASAIR 
WRITE(*,*) "Total water side surface area (sq. in.):" FORMAT (//,"**HEAT EXCHANGER MODEL - OVERALL COEFFICIENT***,//, 
READ *,ASWTR +"Enter the following parameters in the units",/, 
WRITE(*,*) "Fin surface area/total air side surface area:" +"indicated:",/) 
READ *,ASFASAIR FORMAT (" #",I2," Convection coefficient - air side", 
+" (BTU/hr sqft R):") 
WRITE(*,*)"“Enter ‘1’ if unit heater or ‘2’ if duct heater:”" FORMAT (////, 
READ *, FLAG +T25, "HEAT EXCHANGER OVERALL COEFFICIENT’, /, 
IF (FLAG .EQ. 1)THEN +T25, " 
WRITE(*,*) "Conduction coefficient - tube (BTU/hr ft R):” +T5, "Heat exchanger number:",////, 
READ *,KTUBE +T8, "H-WTR = cay BLUsRY Sart RK"; 7, 
WRITE(*,*} "Tube thickness (inches) :" +T8, “KFIN = sh," BIUsShy Et R*:.7, 
READ *, TTUBE +T8, "LFIN = me inches", /, 
WRITE(*,*) "Water side surface area of one tube (sq. in.):" +T8, "WFIN = .3,”" inches",/ 
READ *,ATUBE +T8, “TFIN = af, inches", /, 
WRITE (55,6)HWIR, KFIN, LFIN, WEIN, TFIN, ASAIR, ASWTR, ASFASAIR, +T8, "ASAIR = .1,”" square inches",/, 
+ KTUBE, TTUBE, ATUBE, RFAIR, RFWTR +T8, “ASWTR = Pex square inches",/, 
ENDIF +T8, "ASF/ASAIR = ay Pe ar 
IF(FLAG .EQ. 2)THEN +T8, "KTUBE = .1," BTU/hr ft R",/, 
WRITE (*,*) "Conduction coefficient - tube (BTU/hr ft R):" +T8, "TTUBE = .3," inches", /, 
READ *,KTUBE +T8, "ATUBE = .1,"* square inches", /, 
WRITE(*,*)} "Length of one tube {inches):" +T8, “RFAIR = -4," hr sqft R/BTU",/, 
READ *,LTUBE +T8, "RFWTR er -4,"% hr sqft R/BTU",///, 
WRITE (*,*) "Tube outer diameter (inches):" +T5," H-AIR ETA 1 FIN ETA OVERALL U-AIR",/, 
READ *, ODTUBE +T5, "(BTU/hr sqft R)} 
WRITE(*,*) "Tube inner diameter (inches):" +T5," 
READ *, IDTUBE FORMAT (////, 
WRITE (55, 7)HWTR, KFIN, LFIN, WEIN, TFIN, ASAIR, ASWTR, ASFASAIR, +T25, "HEAT EXCHANGER OVERALL COEFFICIENT", /, 
+ KTUBE,LTUBE, ODTUBE, IDTUBE, RFAIR, RFWTR 
ENDIF +T5,"Heat exchanger number:",////, 
+T8, "H-WTR = ",F6.1," BTU/hr sqft R",/, 
Begin DO loop to compute parameters for each air side convection +T8, "KFIN 2.7 FS 2" BIU/by £E°R"5 7; 
coefficient input. +T8, "LFIN = ",F5.3," inches",/, 
+T8, "WFIN = ",F5.3," inches", / 
DO 22 L = 1, N, i +T8, "TFIN = ",F6.4," inches",/, 
+T8, "ASAIR = ",F7.1," square inches",/, 
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+T8, "ASWTR 


| +T8, "ASF/ASAIR 
| +T8, "KTUBE = BTU/hr ft R",/, 
+T8, "LTUBE : inches", /, 
+T8, “ODTUBE 7 inches", /, 
+T8, "IDTUBE = “ inches",/, 
+T8, "RFAIR “ hr sqft R/BTU",/, 
+T8, "RFWTR " hr sqft R/BTU",///, 
+T5," H-AIR ETA 1 FIN ETA OVERALL U-AIR",/, 
+T5,"“{BTU/hr sqft R) (BTU/hr sqft R)",/, 
+T5," ea | 
8 FORMAT(T10,F4.1,7T23,F4.2,7T37,F4.2,T751,F4.1) 
END 


col 














HEAT EXCHANGER OVERALL COEFFICIENT 


de ee en Se ee 
Set ier Fee ie ee et ee ee ee ee ee ee oe ee oe ee eee ee ee 


Heat exchanger number: B-25 (2-73-1 and 2-80-1) 


H-WTR = 4053.0 BTU/hr sgft R 

KFIN = 26.0 BTU/hr ft R 

LFIN = 0.460 inches 

WFIN = 3.510 inches 

TFIN = 0.018 inches 

ASAIR = 3013.0 square inches 

ASWTR = 465.0 square inches 

ASF/ASAIR = 0.85 

KTUBE = 26.0 BTU/hr ft R 

TTUBE = 0.060 inches 

ATUBE = 67.0 square inches 

RFAIR = 0.0020 hr sqft R/BTU 

RFWTR = 0.0005 hr sqft R/BTU 

H-AIR ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr sqft R) (BTU/hr sqft R) 

6.0 0.87 0.89 4.9 
7.0 0.85 0.87 5.6 
7.9 0.83 0.86 6.1 
8.3 0.83 0.85 6.4 
8.7 0.82 0.85 6.6 
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HEAT EXCHANGER OVERALL COEFFICIENT 


Heat exchanger number: B~-70 (2-8-0, 2-40-1, 2-40-2, 2-60-1, 2-60-2) 


H-WTR 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
TTUBE 
ATUBE 
RFAIR 
RFWTR 


4053.0 BTU/hr sqft R 
26.0 BTU/hr ft R 
0.460 inches 
3.510 inches 
0.0180 inches 

5198.0 square inches 
807.0 square inches 
0.85 

26.0 BTU/hr ft R 
0.060 inches 

90.0 square inches 
0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


fying bd en th nk 


H-AIR ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr sqft R) (BTU/hr sqft R) 


LOI 














Heat exchanger number: 


ASF/ASAIR 


eu haut ttn tana 





HEAT EXCHANGER OVERALL COEFFICIENT 


ed de eee 
Seen ee ee ee ee 


1330.0 BTU/hr sqft R 


01-24-1 


BTU/hr ft R 


inches 
inches 
inches 


17700.0 square inches 
1326.0 square inches 


BTU/hr ft R 


inches 
inches 
inches 


0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


ETA 1 FIN 


(BTU/hr sqft 


ETA OVERALL 
(BTU/hr sqft R) 


le ee ld ee en enn, ee eee 
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HEAT EXCHANGER OVERALL COEFFICIENT 


ee ee et te a a ee ee ee ee ee ewe ee ee ee ee ee ee 
er 


Heat exchanger number: 01-25-1 


H-WTR = 4165.0 BTU/hr sqft R 
KFIN = 118.0 BTU/hr ft R 
LFIN = 0.438 inches 
WFIN = 1.299 inches 
TFIN = 0.0065 inches 
ASAIR = 7959.0 square inches 
ASWTR = 335.0 square inches 
ASF/ASAIR = 0.96 
KTUBE = 221.0 BTU/hr ft R 
LTUBE = 24.0 inches 
ODTUBE = 0.625 inches 
IDTUBE = 0.555 inches 
RFAIR = 0.0020 hr sqft R/BTU 
RFWTR = 0.0005 hr sqft R/BTU 
H-AIR ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr sqft R) (BTU/hr sqft R) 
A nee 
o> 8.2 0.90 0.90 6.4 
WO 8.9 0.89 0.90 6.8 
9.2 0.89 0.89 6.9 
10.90 0.88 0.88 7.4 
10.7 0.87 0.88 7.8 
11.2 0.87 0.87 8.0 





Heat exchanger number: 





HEAT EXCHANGER OVERALL COEFFICIENT 


eC a — er ee i ee ee ee ie oe ss ee ee ee 
a ee 


4165.0 BTU/hr sqft R 
118.0 BTU/hr ft R 


0.438 inches 
1.299 inches 
0.0065 inches 


ASF/ASAIR 


| | | | 


5969.0 square inches 
251.0 square inches 


01-50-0 


BTU/hr ft R 


inches 
inches 
inches 


0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


ETA 1 FIN 


(BTU/hr sqft 


ETA OVERALL 
(BTU/hr sqft R) 


— me ee eee ne ee eee eee ee re em mets ee ee ee ee ee ee ee ee ee se 
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HEAT EXCHANGER OVERALL COEFFICIENT 


Heat exchanger number: 2-25-1 


4165.0 BTU/hr sqft R 
118.0 BTU/hr £t R 
0.438 inches 

1.299 inches 

0.0100 inches 

2868.0 square inches 
230.0 square inches 
0.92 

221.0 BTU/hr ft R 
16.5 inches 

0.625 inches 

0.555 inches 

0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
LTUBE 
ODTUBE 
IDTUBE 
RFAIR 
RFWTR 


it 


ned ok dt wt aon wt ot ot 


H-AIR ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr sqft (BTU/hr 


[TT 








HEAT EXCHANGER OVERALL COEFFICIENT 


ee ee ee 
See eS ee Se eee ee 


Heat exchanger number: 2-16-1 


H-WTR 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 


ASF/ASAIR 


KTUBE 
LTUBE 
ODTUB 


E 


IDTUBE 


RFAIR 
RFWTR 


K-A 
(BTU/hr 


Sl ll ee ee el ee hen hae ee nd 2 


IR 
sqft 


Hu HHH HU eb ot We 


4165.0 BTU/hr sqft R 


118.0 
0.438 
1.299 
0.006 


BTU/hr ft R 
inches 
inches 
5 inches 


3979.0 square inches 
167.0 square inches 


0.96 
221.0 
12.0 
0.625 
0.555 


BTU/hr ft R 
inches 
inches 
inches 


0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


1 FIN ETA OVERALL 


U-AIR 


(BTU/hr 


sqft R) 
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PROGRAM HEAT EXCHANGER OA COEFFICIENT STEAM 


dQ Mgaanan 


REAL HAIR(20),KFIN, LFIN, WFIN, TFIN, ASAIR, ASWTR, ASFASAIR, 
+KTUBE, TTUBE, ATUBE, LTUBE, ODTUBE, IDTUBE, P,AX,M, LC, ETAF, ETAO, 
+UAIR,RFAIR,RFWIR,PI ° 

INTEGER FLAG,N,J,L 
DATA RFAIR, RFWTR, PI/0.002,0.0005,3.14159/ 








Input water and air side convection coefficients and heat exchanger 
dimensional characteristics. Classify computations as either for 
a unit heater or a duct heater. 







NaAaAAN 





OPEN(55, FILE = ‘RESULTS.HXU’ ) 
WRITE (*,4) 
WRITE (*,*) "Quantity of convection coefficients - air side that" 
WRITE(*,*) "will be entered.” 
READ *,N 
pO 21.7 = 1, N; i 

WRITE (*,5)7 

READ *,HAIR(J) 

21 CONTINUE 
WRITE(*,*)} "Conduction coefficient - fin (BTU/hr ft R}:" 
READ *,KFIN 
WRITE(*,*) "Fin length (inches):” 
READ *,LFIN 
WRITE(*,*} "Fin width (inches):" 
READ *,WFIN 
WRITE(*,*)"Fin thickness (inches):” 
READ *,TFIN 
WRITE(*,*) "Total air side surface area (sq. in.):" 
READ *,ASAIR 
WRITE(*,*) "Total water side surface area (sq. in.):" 
READ *,ASWTR 
WRITE(*,*) "Fin surface area/total air side surface area:”" 
READ *,ASFASAIR 































WRITE(*,*}“Enter ’1’ if unit heater or ‘2’ if duct heater:" 
READ *, FLAG 
IF (FLAG .EQ. 1)THEN 
WRITE (*,*) "Conduction coefficient - tube (BTU/hr ft R):" 
READ *,KTUBE 
WRITE(*,*) "Tube thickness (inches) :" 
READ *,TTUBE 
WRITE (*,*) "Water side surface area of one tube (sq. in.):" 
READ *, ATUBE 
WRITE (55,6) KFIN, LFIN, WEIN, TFIN, ASAIR, ASWTR, ASFASATR, 
+ KTUBE, TTUBE, ATUBE, RFAIR, RFWTR 
ENDIF 
IF(FLAG .EQ. 2)THEN 
WRITE (*,*) "Conduction coefficient ~ tube (BTU/hr ft R):" 
READ *,KTUBE 
WRITE(*,*) "Length of one tube (inches):” 
READ *,LTUBE 
WRITE(*,*) "Tube outer diameter (incheés):" 
READ *, ODTUBE 
WRITE(*,*) "Tube inner diameter (inches):" 
READ *, IDTUBE 
WRITE(55,7)KFIN, LFIN, WFIN, TFIN, ASAIR, ASWTR, ASFASAIR, 
+ KTUBE, LTUBE, ODTUBE, IDTUBE, RFAIR, RFWTR 
ENDIF 

























Begin DO loop to compute parameters for each air side convection 
coefficient input. 







DO 22 L = 1, N, 1 





Using relations defined in Section IV-B-2-c, compute single fin 
efficiency and overall efficiency. 


qaaAaAa ANNAN 















efter nt 


as ; ES . SS ea og Soi ees 







((2.0*WFIN)+(2.0*TFIN) )/12.0 
AX = (WFIN*TFIN) /144.0 

M = SQRT( (HAIR(L)*P) / (KFIN*AX) ) 
LC = (LFIN+(TFIN/2.0))/12.0 

ETAF = (TANH(M*LC) ) / (M*LC) 

ETAO = 1.0-ASFASAIR* (1.0-ETAF) 









Using relations defined in Section IV-B-3, compute overall heat 
transfer coefficient. 






aqangan 


IF (FLAG .EQ. 1)THEN 
UAIR = 1.0/ 














+ (((1.0/ (ETAO*HATR(L) * (ASATR/144.0)))+ 
+ (RFAIR/ (ETAO* (ASAIR/144.0)))+4 
+ (RFWTR/ (ASWTR/144.0))+ 
+ | ((TTUBE/12.0) / (KTUBE* (ATUBE/144.0)))) 
+ * (ASAIR/144.0) ) 

ENDIF 

IF (FLAG .EQ. 2)THEN 

UAIR = 1.0/ 

+ (((1.0/ (ETAO*HAIR(L) * (ASAIR/144.0)))+ 
+ (RFAIR/ (ETAO* (ASAIR/144.0)))+4 
+ (RFWTR/ (ASWTR/144.0))+ 
+ (LOG (ODTUBE/IDTUBE) / (2.0*PI* (LTUBE/12.0) *KTUBE) ) ) 
+ * (ASAIR/144.0)) 





ENDIF 






Print results to output file. 






WRITE (55,8)HAIR(L) , ETAF, ETAO, UAIR 






Continue DO loop to completion. 






22 CONTINUE 


a aagqd aan 


4 FORMAT(//,"**HEAT EXCHANGER MODEL - OVERALL COEFFICIENT**",//, 
+"Enter the following parameters in the units",/, 
+"indicated:",/) 

5 FORMAT(" #"*,I2," Convection coefficient - air side", 
+" (BTU/hr sqft R):") 

& FORMAT(////, 
+T25, "HEAT EXCHANGER OVERALL COEFFICIENT", /, 
+T25, “eeommmersrtstsest erases SSS SS SS eee eee" r 4/4 f/f . 
+T5,"Heat exchanger number:",////, 












+T8, "H-STM = VERY LARGE", /, 

+T8, "KFIN = “",F5.1," BTU/hr ft R",/, 

+T8, “LFIN = ",F5.3,° inches", /, 

+T8, "WFIN = “,F5.3," inches", / 

+T8, "TFIN = ",F6.4," inches",/, 

+T8, "ASAIR = ",F7.1," square inches", /, 

+T8, "ASWTR = ",F6.1," square inches", /, 

+T8, "ASF/ASAIR = ",F4.2,/, 

+T8, "KTUBE = ",F5.1," BTU/hr ft R",/, 

+T8, "TTUBE =“, FS.3." inehes” if, 

+T8, "ATUBE = ",F5.1,"° square inches", /, 

+T8, "RFAIR = ",F6.4," hr sqft R/BTU",/, 

+T8, "RFWTR = ",F6.4," hr sqft R/BTU",///, 

+T5," H-AIR ETA 1 FIN ETA OVERALL U-AIR", /, 
+T5," (BTU/hr sqft R) (BTU/hr sqft R)",/, 
#75, "----------------- +--+ + +--+ ee = ++ - ") 






7 FORMAT(////,  - 
+T25, "HEAT EXCHANGER OVERALL COEFFICIENT", /, 
+725, “saesssssasessrsssssssesesnzeessre=",////, 
+T5,"Heat exchanger number:",////, 







+T8, "H-STM = VERY LARGE", /, 
+T8, "KFIN = ",F5.1,"° BTU/hr ft R",/, 
+T8, "LFIN = ",.F5.3," inches", /, 

+T8, "WFIN = ",FS.3," inches", / 

+T8, "TFIN = ",F6.4," inches",/, 

+T8, "ASAIR = ",F7.1," square inches", /, 
+T8, "ASWTR = 1," square inches",/, 





+T8, “ASF/ASAIR 






= ",F5.1," BTU/hr ft R",/, 

= *,F5.1," inches", /, 
+T8, “ODTUBE = ",F5.3,"“ inches’,/, 
+T§8, “IDTUBE = ",F5.3," inches",/, 
+T8,"RFAIR = ",F6.4," hr sqft R/BTU",/, 
+T8, "RFWTR = ",F6.4," hr sqft R/BTU",///, 


+75," H-AIR ETA FIN ETA OVERALL 
+75, "(BTU/hr sqft R) 


8 FORMAT(T10,F4.1,T23,F4.2,T737,F4.2,T51,F4,1) 
END 


U-AIR",/, 


(BTU/hr sqft R)",/, 
475 , Mem nr ee ee ee ee ee ee 
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HEAT EXCHANGER OVERALL COEFFICIENT 


Heat exchanger number: B-25 (2-73-1 and 2-80-1) 


VERY LARGE 

26.0 BTU/hr ft R 
0.460 inches 
3.510 inches 

0.0180 inches 

3013.0 square inches 
465.0 square inches 
0.85 

26.0 BTU/hr ft R 
0.060 inches 

67.0 square inches 
0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


H-STM 
KFIN 
LFIN 
WIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
TTUBE 
ATUBE 
RFAIR 
RFWTR 


pnenk t &nrnrnt tous 


H-AIR ETA 1 FIN ETA OVERALL U~-AIR 
(BTU/hr sqft (BTU/hr sqft R) 


Sil 





Heat exchanger 


H-STM 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF /ASAIR 
KTUBE 
TTUBE 
ATUBE 
RFAIR 
RFWTR 


H-AIR 
(BTU/hr sqft 


0.85 


peuu dt? iungne tb ut 


ETA 1 FIN ETA OVERALL 
R) 
0.82 0.84 
0.81 0.84 
0.80 0.83 
0.80 0.83 
0.79 0.83 
0.79 0.82 





HEAT EXCHANGER OVERALL COEFFICIENT 


oe oe 


ee ey ie re ee ee ee ee ee ee ee 


number: 


B-70 (2-8-0, 2-40-1, 2-40-2, 2-60-1, 2-60~2) 


VERY LARGE 
26.0 BTU/hr ft R 
0.460 inches 
3.510 inches 
0.0180 inches 
5198.0 square inches 
807.0 square inches 


26.0 BTU/hr ft R 
0.060 inches 

90.0 square inches 
0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


U-AIR 
(BTU/hr 


et ee cee eee me ce ee ee rere wre eet ee mee ee eee eee eee ee ee ee 


sqft R) 
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Heat exchanger 


H-STM 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
LTUBE 
ODTUBE 
IDTUBE 
RFAIR 
RFWTR 


H-AIR 
(BTU/hr sqft 


Vn n t tot dt wt ow ke & nn 





HEAT EXCHANGER OVERALL COEFFICIENT 


number: 01~24-1 


VERY LARGE 

118.0 BTU/hr ft R 
0.450 inches 

1.880 inches 

0.0080 inches 
17700.0 square inches 
1326.0 square inches 
0.91 

221.0 BTU/hr ft R 
34.0 inches 

1.125 inches 

1.035 inches 

0.0020 hr sqft R/BTU 
0.0005 hr sqft R/BTU 


ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr 

















HEAT EXCHANGER OVERALL COEFFICIENT 


ee Baer eed a ee 
we ee Oe es Oe et ee ee oe a a 


Heat exchanger number: 01-25-i1 


H-STM = VERY LARGE 
KFIN = 118.0 BTU/hr ft R 
LFIN = 0.438 inches 
WFIN = 1.299 inches 
TFIN = 0.0065 inches 
ASAIR = 7959.0 square inches 
ASWTR = 335.0 square inches 
ASF/ASAIR = 0.96 
KTUBE = 221.0 BTU/hr ft R 
LTUBE = 24.0 inches 
ODTUBE = 0.625 inches 
IDTUBE = 0.555 inches 
RFAIR = 0.0020 hr sqft R/BTU 
RFWTR = 0.0005 hr sqft R/BTU 
H-AIR ETA i FIN ETA OVERALL U-AIR 
(BTU/hr sqft R) (BTU/hr sqft R) 
8.2 0.90 0.90 6.6 
8.9 0.89 0.89 7.0 
9.2 0.89 0.89 wae 
100 0.88 0.88 i Pe 
10.7 0.87 0.88 8.1 
LL.2 0.87 0.87 8.4 
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HEAT EXCHANGER OVERALL COEFFICIENT 


Heat exchanger number: 01-50-0 


VERY LARGE 

118.0 BTU/hr ft R 
0.438 inches 

1.299 inches 

0.0065 inches 
5969.0 square inches 
251.0 square inches 

0.96 

221.0 BTU/hr ft R 
18.0 inches 

0.625 inches 

0.555 inches 

0.0020 hr sqft R/BTU 

0.0005 hr sqft R/BTU 


H-STM 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
LTUBE 
ODTUBE 
IDTUBE 
RFAIR 
RFWTR 


tuerkt Font ir tnt 


H-AIR ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr saft (BTU/hr sqft R} 


6lt 

















HEAT EXCHANGER OVERALL COEFFICIENT 


SS rm ee ee ee ee a i i a en ee ae a a ee 
ee a a a ee a ee eS 


Heat exchanger number: 2-25-1 


H-STM = VERY LARGE 

KFIN = 118.0 BTU/hr ft R 

LFIN = 0.438 inches 

WFIN = 1.299 inches 

TFIN = 0.0100 inches 

ASAIR = 2868.0 square inches 

ASWTR = 230.0 square inches 

ASF/ASAIR = 0.92 

KTUBE = 221.0 BTU/hr ft R 

LTUBE = 16.5 inches 

ODTUBE = 0.625 inches 

IDTUBE = 0.555 inches 

RFAIR = 0.0020 hr sqft R/BTU 

RFWTR = 0.0005 hr sqft R/BTU 
H-AIR ETA 1 FIN ETA OVERALL U-AIR 

(BTU/hr sqft R) (BTU/hr sqft R) 

10.1 0.92 0.92 8.5 
1 ew 0.91 0.92 9.3 
12.1 0.90 0.91 9.9 
12.9 0.89 0.90 10.5 
L338 0.89 0.90 y iy a 8 
14.5 0.88 0.89 11.5 
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HEAT EXCHANGER OVERALL COEFFICIENT 


Heat exchanger number: 2-16-1 


H-STM 
KFIN 
LFIN 
WFIN 
TFIN 
ASAIR 
ASWTR 
ASF/ASAIR 
KTUBE 
LTUBE 
ODTUBE 
IDTUBE 
RFAIR 
RFWTR 


H-AIR 
(BTU/hr 


va 


nt nbwtkt t t bt ne vob 


VERY LARGE 

118.0 BTU/hr ft R 
0.438 inches 

1.299 inches 

0.0065 inches 
3979.0 square inches 
167.0 square inches 
0.96 

221.0 BTU/hr ft R 
12.0 inches 

0.625 inches 

0.555 inches 

0.0020 hr sqft R/BTU 

0.0005 hr sqft R/BTU 


ETA 1 FIN ETA OVERALL U-AIR 
(BTU/hr sqft R) 





Cc 
Cc 
Cc 
Cc 
c 
Cc 


aqanNaAn 


qanqgaaQg Aa 


AaAaAQ 


ANNAN 
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PROGRAM HEAT EXCHANGER MODEL NTU ANALYSIS 


REAL T3,T1,VDOTAIR, UAIR, ASAIR, QDOT, MDOTAIR, DAIR, T, PATM,R, 


+CAIR, CPAIR, MDOTWTRI, MDOTWTR, CWTR, CPWTR, QDOTMAX, CMIN, CMAX, CR, 
+EFF,NTUIL, EFFI,NTU,ASAIRI, VDOTWTR, T2,T4, QDOTFINAL 


INTEGER FLAG 
DATA CPAIR, CPWTR, TSTD/0.24,1.0025,70.0/ 


Input water and air inlet temperatures, overall coefficient, air flow 
rate, air side surface area, and heat transfer rate. 


OPEN(55, FILE = ’RESULTS.HXNTU’ ) 


WRITE (*, 4) 

WRITE(*,*) "Water entering temperature (degrees F):" 

READ *,T3 

WRITE(*,*)"“Air stream entering temperature (degrees F):" 
READ *,T1 


WRITE(*,*) “Air stream volumetric flow rate (CFM):" 

READ *,VDOTAIR 

WRITE (*,*) "Overall heat exchanger coefficient (BTU/hr sqft R):" 
READ *,UAIR 

WRITE(*,*} "Total air side surface area (sq. in.):" 

READ *,ASAIR 

WRITE(*,*) "Heat transfer rate required (BTU/hr):" 

READ *,QDOT 


Compute air mass flow rate and heat capacity rate. 


MDOTAIR = DAIR(TSTD) *VDOTAIR*60.0 
CAIR = MDOTAIR*CPAIR 
WRITE (55,5)T3,T1,VDOTAIR, UAIR, ASAIR, QDOT, MDOTAIR, CAIR 


MDOTWTRI = QDOT/ (CPWTR* (T3-T1) } 


Begin DO loop to compute Effectiveness-NTU parameters at water mass 
fiow rates iterated up to 4876 lbm/hr (10 gpm). 


DO 23 MDOTWTR = MDOTWTRI, 4876.0, 10.0 
IF(MDOTWTR .GT. 4866.0) THEN 
WRITE (55,8) 
GO TO 99 
ENDIF 


Compute water heat capacity rate, CWTR, and compare with air heat capacity 
rate, CAIR. If CAIR is less than CWTR, proceed with CMIN equal to CAIR. 


CWIR = MDOTWTR*CPWTR 
IF(CAIR .LT. CWIR) THEN 
CMIN CAIR 
CMAX CWTR 
FLAG ul 
QDOTMAX = CMIN* (T3-Ti} 
CR = CMIN/CMAX 
EFF = QDOT/QDOTMAX 
DO 20 NTUI = 0.0, 5.5, 0.01 
EFFI = 1.0 - EXP((1.0/CR)* (NTUI**0.22) * 


toast 


+ (EXP (-CR* (NTUI**0.78))-1.0)) 


IF (ABS (EFFI-EFF) .LT. 0.01) THEN 
NTU = NTUI 

T2 = Tl + (QDOT/CAIR) 

T4 = T3 - (QDOT/CWTR) 


After computations of heat capacity ratio, CR, effectiveness, EFF, 
NTU, T2, and T4, proceed to line 22. 


GO TO 22 
ENDIF 
CONTINUE 


C 
C 
Cc 


A000 


ANAAN 


9 ANAAN 


21 


22 


23 
4 


5 


6 





If CWTR is less than CAIR, proceed with CMIN equal to CWTR. 


IF(CWTR .LT. CAIR) THEN 
CWTR 


CMIN = 
CMAX = CAIR 
FLAG = 2 


QDOTMAX = CMIN* (T3-T1) 
CR = CMIN/CMAX 
EFF = QDOT/QDOTMAX 
DO 21 NTUI = 0.0, 5.5, 0.01 
EFFI = 1.0 - EXP((1.0/CR) * (NTUI**0.22)* 


+ (EXP (-CR* (NTUI**0.78))-1.0)) 
IF(ABS(EFFI-EFF) .LT. 0.01)THEN 
NTU = NTUI 
T2 Tl + (QDOT/CAIR} 


T4 T3 ~ (QDOT/CWTR) 
After computations of heat capacity ratio, CR, effectiveness, EFF, 
NTU, T2, and T4, proceed to line 22. 


GO TO 22 


Compute heat exchanger air side surface area and compare with actual 
air side surface area. If computed value is less than actual value, 
print results to output file. 


IF(NTU .LT. 1.0E~-20)NTU=5.0 
ASAIRI = ((NTU*CMIN) /UAIR) *144.0 
VDOTWTR = (MDOTWTR*7.481)/({(60.8*60.0) 
QDOTFINAL = EFF*CMIN* (T3-T1) 
WRITE (55,6)MDOTWTR, VDOTWTR, CMIN, CMAX, CR, EFF, NTU, ASAIRI 
IF(ASAIRI .LT. ASAIR) THEN 

WRITE (55, 7)MDOTWTR, VDOTWTR, CMIN, CMAX, CR, EFF, NTU, ASAIRI, 

+ T2,T4, QDOTFINAL 
IF(FLAG .EQ. 2) THEN 
MDOTWTRI = CAIR/CPWTR + 1.0 


GO TO 19 
ENDIF 
IF(FLAG .EQ. 1)GO TO 99 
ENDIF 


If computed air side surface area is greater than actual, continue 
DO loop until solution is found or end program at water mass flow rate 
of 4,876 lbm/hr (10 gpm). 


CONTINUE 
FORMAT (//,"****HEAT EXCHANGER MODEL - NTU ANALYSIS****",//, 


+"Enter the following parameters in the units", /, 
+"indicated:",/) 


FORMAT (////, 
+T25, "HEAT EXCHANGER MODEL - NTU ANALYSIS",/, 
4725, “eesseressssssassssscsrccessssesese==",////, 


+T5, "Heat exchanger number:",////, 


+T8,"T3 = ",F5.1," degrees F",/, 

+78, "T1 = ",FS5.1," degrees F",/, 

+T8,""VDOTAIR = ",F6.1," CFM", /, 

+T8, “UAIR = ",F5S,1,;" BTU/hr sqft. R";/; 

+T8, "ASAIR = ",F7.1," square inches", /, 

+T8, "“QDOT = ",F8.1," BrU/hr",/, 

+T8, "MDOTAIR = ",F7.1," ibm/hr’, /, 

+T8, "CAIR me 96 CL,” BTU/hr R .7//S/, 

+T1," Mass Volume Cc Cc CR Ef£ NTU Hx Surface",/, 
+T1," Flow Flow min max Area’,/, 
+T1,"(lbm/hr) (gpm) (BTU/hr R) (sq. in.)",/, 
+T1, OU is Facets cclas Ste Sg eet Sets a er “as me Ps ss Sa cana ec cack eel ems a at es "ea ccs ls esa ces ls cn susan Tes “ } 


FORMAT (/,T2,F6.0,T10,F5.2,T16,F6.1,T724,F7.1,T33,F4.2,T39,F4.2, 
+T45,F4.2,T52,F6.0) 


199 


ecl 































+T1,". Mass Volume ie Cc CR Eff NTU Hx Surface", /, 


+T1," Flow Flow min max Area",/, 
+T1,"{lbm/hr} (gpm) (BTU/hr R) (sq. in.)",/, 
HTL, "~~-------------- =o en nnn nnn ne nee eee eee = "I; 


+7T2,F6.0,T10,F4.2,T16,F6.1,T724,F7.1,T33,F4.2,T739,F4.2, 
+T45,F4.2,T52,F6.0,//, 

+T3,°T2 ",F5.1," degrees F",/, 

+T3, "T4 ",F5.1," degrees F",/, 

+T3, "“QDOTFINAL " F8.1," BTU/hr", //, 
+////1) 
8 FORMAT(//,T5,"At water flow rates of up to 100 gpm,”,/, 
+T5, "there is no solution found.") 









ee tt dt. 
zs SSS SSS SS SS SS SSS SSS ee Oe eee 






FUNCTION subroutines to compute temperature dependent transport 
property of air. 


AaAnNAAA 


FUNCTION DAIR(T) 

REAL DAIR,T,PATM,R 

DATA PATM,R/14.7,53.34/ 

DAIR = (PATM*144.0)/(R* (T+459.67)) 
RETURN 

END 





HEAT EXCHANGER MODEL - NTU OPTIMIZATION ANALYSIS 


ee 


Heat exchanger number: 


me re me cp cm em ek ates ee meme my ese ee meee eet i te eee ee ee ee ee ee ee 


T3 
T1 
VDOTAIR 
UAIR 
ASAIR 
QDOT 
MDOTAIR 
CAIR 
Mass Volume 
Flow Flow 
(lbm/hr) (gpm) 
402. 0.82 
412. 0.85 
422. 0.87 
432. 0.89 
442. 0.91 
452. 0.93 
462. 0.95 
472. 0.97 
482. 0.99 
492, 1.01 
502. 1.03 
S12. 1.05 
522. 1.07 
532. 1.09 
542. a oe De @ 
552. 1.13 
562. 1.15 
572. 1,i7 
582. 1.19 
592. 1.21 


vuni t t a & t 


190.0 degrees F 
62.0 degrees F 


3680.0 CFM 


4.8 BTU/hr sqft R 
17700.0 square inches 
51600.0 BTU/hr 
16543.2 lbm/hr 


01-24-1 


3970.4 BTU/hr R 


Cc Cc 
min max 
(BTU/hr R) 
403.1 3970 
413.2 3970 
423.2 3970 
433.2 3970 
443.2 3970 
453.2 3970 
463.3 3970 
473.3 3970 
483.3 3970 
493.4 3970 
503.4 3970 
513.4 3970 
523.4 3970 
53330 3970 
543.5 3970 
J33:5 3970 
563: 5 3970 
573. 3370 
583.6 3970 
593.6 3970 


r FF Fe FF F&F FB LH fF 


> fb bk f 8 fF fk ek BP fe be Oe 


a] 


Eff NTU Hx Surface 


oe — Sere ee 
a Ic a 


Area 
(sq. in.) 


124 


632. 1.30 633.7 3970.4 16 0.64 1.07 20342 
642. 1.32 643.7 3970.4 16 0.63 1.05 20277 
652. Les -653.7 3970.4 16 0.62 1.02 20005 
662. 1.36 663.8 3970.4 0.17 0.61 0.99 19714 
672. 1.38 673.8 3970.4 0.17 0.60 0.97 19608 
682. 1.40 683.8 3970.4 0.17 0.59 0.94 19284 
692. 1.42 693.9 3970.4 0.17 0.58 0.92 19150 
702. 1.44 703.9 3970.4 0.18 0.57 0.90 19005 
7125 1.46 713.9 3970.4 0.18 0.56 0.88 18847 
722 1.48 723.9 3970.4 0.18 0.56 0.86 18677 
732 1.50 734.0 3970.4 0.18 0.55 0.84 18496 
742. 1.52 744.0 3970.4 0.19 0.54 0.83 18525 
Tes 1.54 754.0 3970.4 0.19 0.53 0.81 18322 
762 1.56 764.0 3970.4 0.19 0.53 0.79 18107 
772 1.58 774.1 3970.4 0.19 0.52 0.78 18113 
782 1.60 784.1 3970.4 0.20 0.51 0.76 17877 
792. 1.62 794.1 3970.4 0.20 0.51 0.75 17867 
802. 1.64 804.1 3970.4 0.20 0.50 0.74 17852 
812. 1.67 814.2 3970.4 0.21 0.50 0.72 17586. 
f1!!{SOLUTION!!!!!SOLUTION! !!!!!SOLUTION!!!!!SOLUTION! !!! 
Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow min max Area 


(lbm/hr) (gpm) (BTU/hr R) (sq. in.) 


-—— ee ee ee eee ee ee eee ee eee ee ee ee ee ee ee es er eae ee ee 


812. 1.67 814.2 3970.4. 0.21. 0:56 0.72 17586. 


= 75.0 degrees F 
T4 = 126.6 degrees F 
= 51600.0 BTU/hr 


3961. 8.12 3970.4 gofes@ 1.00) “GCG. 10 207.14 13102. 


Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow min max Area 
(lbm/hr) (gpm) {BTU/hr R} (sq. in.) 


——e me te a ee ee ee ee ae ce crt ce ee re ee ee te ee ee eee yee ee ee ee ee 
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3961. 8.12 3970.4 3971.4 1.00 06.10 9.1 


T2 = 75.0 degrees F 
T4 177.0 degrees F 
QDOTFINAL = 51600.0 BTU/hr 





HEAT EXCHANGER MODEL - NTU OPTIMIZATION ANALYSIS 


——— 
I - 


Heat exchanger number: 2”°&-o 


T3 = 190.0 degrees F 

Ti = 60.0 degrees F 

VDOTAIR = 1505.0 CFM 

UAIR = 3.7 BTU/hr sqft R 

ASAIR = 5198.0 square inches 

QDOT | 5294.0 BTU/hr 

MDOTAIR = 6765.6 lbm/hr 

CAIR = 1623.8 BTU/hr R 
Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow min max Area 

(lbm/hr) (gpm) (BTU/hr R) (sq. in.) 


eee ce ee me me me re ee te ee ree ee ee ee ee ee ee ee eye ee ee ee oe 


41. 0.08 40.7 1623.8 0.03 1.00 4.81 7623. 
51. 0.10 90.7 1623.8 0.03 0.80 1.61 3180. 


II l'IiSoLUPTION!!!!t !SOLUTION!!!1)}! '!SOLUTION! !! ! '!SOLUTION! !!: 
Mass Volume Cc c CR Eff NTU Hx Surface 
Flow Flow min max Area 

(lbm/hr) (gpm) (BTU/hr R) (sq. in.) 


ee ee oe ee ee ee ee ee eee eee Bee ee ee ew ee ee wee ew ewe ee ee ee ee ee 


T2 = 63.3 degrees F 
T4 = 85.7 degrees F 
QDOTFINAL = 5294.0 BTU/hr 


PLi!!SOLUTION! !!!!SOLUTION! !!!!!SOLUTION!!!!!SOLUTION!!!! 


1621. 3.32 1623.8 1624.8 1.00 0.03 9.02 1264. 


ill! L SOLUTION!!! ! !SOLUTION! $!5§!!SOLUTION! !!! !SOLUTION! !!! 
Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow Min max Area 

{lbm/hr) (gpm) (BTU/hr R) (sq. in.) 


me a mm wee ee ee ree ee ree me mer ree te ee wre ret tee eee eee ee ee eee ee ee ee ee 


1621. 3.32 1623.8 1624.8 1.00 0.03 0.02 1264. 


T2 = 63.3 degrees F 
T4 = 186.7 degrees F 
QDOTFINAL = 5294.0 BTU/hr 





hers 


HEAT EXCHANGER MODEL - NTU OPTIMIZATION ANALYSIS 


Se eed 
=.= SSS SS SS SS SS SS Se eS SS Se SS Se Se ee ee ee ee re ee oe 


Heat exchanger number: 2-80-1 


T3 = 190.0 degrees F 

T1 = 60.0 degrees F 

VDOTAIR = 400.0 CFM 

UAIR = 4.2 BTU/hr sqft R 

ASAIR = 3013.0 square inches 

QDOT = 7747.0 BTU/hr 

MDOTAIR = 1798.2 lbm/hr 

CAIR = 431.6 BTU/hr R 
Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow min max Area 

(lbm/hr} (gpm) (BTU/hr R) (sq. in.) 


ee ee ee eee eee ee Se See eS ee eee Eee eS eee = 


59. 0.12 59:,6 431.6 0.14 1.00 5.00 10216. 
69. 0.14 69.6 431.6 0.16 0.86 2.16 5156. 
79, 0.16 79.6 €31.6: /0.18. ‘G.:75': “1.52 4150. 
89. 0.18 89.7 431.6 0.21 0.66 1.20 3689. 
99: 0.20 ee 431.6 0.23 0.60 1.00 3418. 
109. O.22 109.7 431.6 0.25 0.54 0.86 3235. 
119. 0.24 119.7 431.6 0.28 0.50 0.75 3079. 


129s 0.27 129.8 431.6 0.30 0.46 0.67 2981. 


tIttI SOLUTION!!! ! 'SOLUTION! !!!!!SOLUTION! !!!!SOLUTION! !!! 
Mass Volume Cc Cc CR Eff NTU Hx Surface 
Flow Flow min max Area 

(lbm/hr) (gpm) (BTU/hr R) (sq. in.) 


ee i i ee einai 


129. 0.27 129.8 431.6 0.30 0.46 0.67 2981. 
78.0 degrees F 


T4 = 130.3 degrees F 
= 7747.0 BTU/hr - 


431. 0.88 431.6 432.6 1.00 0.14 90.16 2367. 








Mass Volume Cc Cc 
Flow Flow min max 
(lbm/hr) (gpm) (BTU/hr R) 


431.6 432.6 1.00 0.14 0.16 


72 78.0 degrees F 
T4 172.1 degrees F 
QDOTFINAL 7747.0 BTU/hr 
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APPENDIX D. EFFECTIVENESS - NTU ANALYSIS 
TABULAR RESULTS 
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Ocl 


> 
oOo = 
ro) 
ae 
- 
= 
bl 


187 
187 


1 CO 
1 0 
- 


79 


16270 
16722 


aaa 
Fi Saal 
ey! 
tea, 
[ed 
a 
eel 
fae ase! 
ae 
a 
stated 
| 12 

87 

85 

| 14807 | 83 | 187 | 


CR 


omen 
= 
a 
Reece tees 
= 
CR 
0.21 
0.19 
| 0.13 | 048 | 





err) 
Wee al 
ea 
He atl 
i 
Ee 
ow 
ae 
NTU 


QMAX 





CWTR 


| 4886 | 98149 | O18 | 015 | 016 | 15663 [| 81 | 187 | 
0.18 


Peden 28 
WATER PROPERTIES 
T3 
190 
fe aol 
ae 
4886 
eee 
Pode 
oe! 
ey 
Paar 
sae 
CWIR 


CP 
74008 
a 

4876 
oes 
a 
3 
alate 
a 





Rises toed 

| 

= 
T1 

L ie0e = 

ae 

eet 
CAIR 
431 

| 539 | 
647 
755 

Pie est 
CAIR 

















K 
0.0148 


MU 
1.23E-05 
ir 
ae 
H-AIR 
| 6.0 
7.0 
7.9 
8.7 
Pegecsel 
| 
aioe 
a 
ead 
oe 
ees 
H-AiR 
ene 
AiR PROPERTIES 
MU 
1.23E-05 
are 
gl tio 
H-AIR H-STM 
>>Q 


| 
Pee os esl 
AIR PROPERTIES 
CP 
es 
JH 
| 0.0042 | 
aaa 
ae 
a 
JH 
0.0048 
CP 
0.24 
JH 
0.0048 


5962 


= 
bee 
zw 
wv 
: 


Ys) slo = 
Fis] | SIS mi 
= 


roe el 
Poel 
eee 
ee el 
iia 
Mean 2 
Ral 
RHO 
0.0749 
peo ll 
Paes | 
0.0749 
el 
Heoeeee eal 
RE 
2981 





fad 
rr 
rs 
oe 
—, 
DH 
fe — ae 
Hee 
G 
1.16 
2.03 
2.32 
ae 
ee 
= 





MDOTAIR 
1798 


73013 | 


COIL DIMENSIONS 








UNIT HEATER B-25 (2-73-1, 2-80-1 
a 
eo 
COIL DIMENSIONS 
| 600 | 2696 | 
ae 
ne! 
an 
HYDRONIC WITH VARIATION OF T3 
so 
— 4 
nee 
aie 








HYDRONIC 


CFM 
700 
800 

amma 
a 








Tel 
















































































































EO re) Gt OEE A 
aa NPR a! a 
nf 
Sele a Sa 
ree pe ae 
eo eer) ae baa 
a aR Resa esa ieee 
a See ee 
| AFF | cP. | 73 | ae ii oan Paes 
| 120 0.24 | 1236-06 ——— fee 
—e aun rena eee ae NE aes eens 
=a ee eae Ie eee a page te 
CFM G | RE | JH | HAIR CAIR QM NTU Eps | Q | 72 14 
ease | as sor | a 6.7 | 1623 — ae 4886 0.15 | 033 | 0.13 
270 | "0.0033 98.6 _| ‘0 | 1941 043 | 040 | 0.12 
P2100 _| 3.15 ; ai a ee oe ee 
360) | 9241 | 0.0026 | 10.1 | 4053 | 7. 4886 | 336511 | 0.10 | 053 | 0.09 | 30834 | 72 | 184 _| 
2700 ) [10396 ~0.0024 | 10. . 378675 | 0.09 | 060 | 0.08 | 31787 | 71 | 184 
4053 | 7.5 | 9236 ee ee Oe ee 
a a ie ieee Bae 
ay +} 
— anel SEE Sana mae! ears 
HYDRONICWITHVARIATIONOFT3 | ia aa eee 
a Sree Se aCe Hoses ches omer AN 
Bl | =  —  — 
Le eet ease | 
pe Ee: ieee Ce See 
. ae ee are aa 
ie ee ae | CR | EPS | 
225 | | —879s 4053 | 6.7 0.33 | 0.13 | 31469 
mat ee a a Seen ieee Sea 
ee en 
a 
es ene: rea eta 
COIL DIMENSIONS ame 
—20__S196_| 088 |__| arte |_O24 
—- — 
ae as: ost aes 
CFM ——.. BE. AR Borst CSTM CR | Eps | Q | BB | 
| 1505 | 6763_ | 2.25 | | 5795 | 0.00387 | 9.0 | | 42.5 | INF. | 292182 _ 0.00 | 0.14 | 40445 | 85 | 



















































































































































DUCT HEATER 01-24-1 ee eS See eee [ —— ————————E 
Sean ines Gaaaor ame sseunt arc a a — oe oes — 
Se —4 a 
aaa Seer MERE aa 
ee eee onsen a Mig Mt 
— — — ao 
| | STAR PROPERTIES _ [WATER PROPERTIES PF 
| AFF | ASAIR | DH [| | RHO ee we 13 - eee ees | 
| 405 | 17700 | 0.25 | | 0.0749 | 0.24 | 1.23E-05/ 0.0148 | 0.72 BZ .. | 4.002 | 190 | es 7 ae 
SS A ON Se SD SRA GORE See: en Soe oe — 
at 
aM MB OTAR 8 —————|-- —| | BE | OAR CAIR Q QMAX | NTU EPS | Q q2 14 
16538 | 163 | | 2762 | «0.004 | 7.0 | 4053 | 57 §08045 | 0.18 | 081 | 0.15 | 74778 81 175 
4500 20223 2.00 3378 0.0036 6.2 4876 4886 | 621251 | 016 | 099 | 0.13 79 
5400 eye ee ee 8.5 5824 4886 0.14 a 172 
6300 e+} ee 9.0 6795 _ 4876 171 
| 7200 | 32357 | 3.20 | | +5404 | 0.0028 | 96 | 4053 | 7.4 | 7766 4876 | 625376 | 0.19 | 063 | 016 | 98424 | 75 | 170 
7880 | $8076 1 _8.2/ | ____1_Se28_1 9.00276 |_a7 _|_4088_{_f4 7938 ; 4876 0.19 | 062 | 016 | 98586 
HYDRONIC WITH VARIATION OF T3 El boo) P| Hes ee nef ee | 
Ser a (ae eee pom eee a 
5 ee pe ee 
ore eeeen Pa 
aa De cet seam lee ea a . Pe ee rn eee —-——~ 
Fee wet ga a | ga “gu awe} geo 
CFM ‘MDOTAIR | JH H-AIR | H-WTR | UAIR MDOTWT | CWTR CR ae > ae ea ae a 7 
| 3680 | 16538 | 1.63 | 2162 0.004 1 TO | A088 | _ $7 | 3980 ____| 4876 | 4886 | 567581 | 018 | 081 | 015 | 83541 | 83 | 188 | 
—— Se poe ae er ae ee eh ee I ete a tle 
a a ae! are sce SGA ANC? SOLS GOLLESE CSGGRET CARAT ORONO ROORAR S IA 
a — | ones 
STEAM ae Eee aan See eg el ee lee od 
| COILDIMENSIONS | | ae = 4 
405 | 17700 | 0.25 |  ~—s|_0.0749 | 0.24 | 1.23E-05] 0.0148 | Oe INF 240 . poe aaa 
Sa a cae Se Sees meee oe oes ae 
ae ike —, -+—— ae Ra seal 
CFM |MDOTAIR| G | | RE | JH | H-AIR UAIR | CR | EPS | T2 
| 3680 | 16538 | 163 |  _| 2762 | 0.0048 | 70 | >>0 [| 61 | 3969 019 | 000 | 0.17 [121615 | 93 | | 
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se 
eee ee el ae ena iat ees Etaen e 
tee ae es, ere, =e — eae, eee 
= ee ee 
i — —_— 
EOL ERUNSONS poof thik brorertes | [aes ae eae SEE eis AIT Ean: mara?! 
[COILDIMENSIONS | | AIRPROPERTIES | | ee ee 
| AFF | ASAIR | DH | | RHO | cp | MU | KO | PR | 1 GR TS a) eaaee! a ree eee 
a ee 
ae FS ase aaa Sai nes meee eres ieee eee! 
Sar eget ge +e 
oe ee | JH | HAIR | HWIR | UAIR | CAIR | MDOTWT| CWIR | QMAX | NTU | CR | EPS | Q | 7 
6292 | 1.62 | 1691 | 0.0047 | 82 | 4053 | 64 | 1510 | | 4876 | 4886 | 332196 | 023 | 0.31 | 0.20 | 66303 | 14 
7640 | 1.97 2054 oe ee | eee te 4876 | 4886 | 403381 | 0.20 | 0.38 | 018 | 71200 | 9 
a ee 0.0037 | 9.2 | 4053 | 69 IE 4876 | 4886 | 474566 | 0.18 | 0.44 73016 | 4 
| 2300 | 10336 | 267 | 2779 | 0.0035 | 10.0 | 4053 | 74 | 2481 4876 > eee | See ee 061 | 0.14 | 78375 | 2 | 174 | 
| 2600 | 11684 | 302 [| 0.0033 7.8 616936 | 0.15 : 0.13 A 173 
j 7800 _| 12888 828 | {0.0032 | 41.2 | 4053 _ | __ 3020 _ 0.15 0.13 
ee ee oe ee a 
eae) eee en ae! ie eee Roe —— 
Pmecemnawalen ce = 
HYDRONIC WITH VARIATION OF Koen — an 
aa ees, eee ae pemmeees| eee ee Sees eee 
———— = pp 
—— pee! 
—- [— —— —- 
Rate 
UAIR | CR | EPS | Q | 72 | 1% | 
| 1400 | 6202 | 1.62 | 6.4 17 190 
ial oon eee pn 
oe ee ae ee ee =p 
a = 
ra ane eeray = eo eel 
| COILDIMENSIONS | | aaraiess, hess Es =! 
AFF ASAIR DH RHO PR T 
a2 pee i PRES 
fe aed od ii 
aie eer ee ee See pane a] tee ea Retin eects 
Gof | RE JH | CHAIR | HSTM | UAIR | CAIR | ~~ MDOTST | CSTM | QMAX | NTU | mz} 
[| 1400 | 6292 | 162 | | 1691 | 0.0047 | 82 [ >>0 | 66 | 1510 [| | 919 | INF. | 407696 | 0.24 | 0.00 | 87408 | 28 | 

































































































DUCT HEATER 01-50-0 — a 
sce ne een Sees ns Sane SN ey Se rae ees a ae ee — 
fiona a Fees anes Lace eee 
HYDRONIC —--—— 
SOIL DIMENSIONS | |_| arn propeRies =P eee 
COILDIMENSIONS || | AIR PROPERTIES _| aan ee an 
| AFF | ASAIR| DH | =| RHO | CP | MU | K | PR | T | hf] cP OU] U3 UT = Mee creel uses 
| 116 | 6969 | 0154 | | 0.0749 | 0.24 | 1.23E-05) 0.0148 | 0.72 | -30 | ~~ | 1.002 | 190 eee sees ee me es 
a Pierre ee ye ee eee Were eben ot (ee Gee (ei eerie 
| CFM |MDOTAIR, G {| | RE | JH | HAIR | HWIR NTU | CR {| EPS | Q | 2 | 74 
a a 0.0056 0 | soo ||: «C4876 | «4886 | 177962 | «2029 | O17 | 024 | 43215 | 23 181 
| 900 | 4045 | 1.39 | | 1453 | 0.0051 | 971 | 4876 | 4886 | 213555 | 0.26 | 0.20 0.22 | 47509 19 180 
| 1050 | 4719 | 163 | | 1695 | 0.0047 | He |. ee 0.23 0.23 | 0.20 | 50341 14 | 180 
1200 5393 1.86 | = | 1937 | 0.0043 | 1294 4876 4886 0.21 0.26 0.18 10 
ae ee eae | 1456 | | 4876 | 4886 | 320332 | 0.19 0.17 | 54312 | 7 | 179 
6741 | 232 | —s| 2421 | 0.0038 | 95 | eis |_|_ 4878 4886 jsssus [oe ons | one [meow | fre 
Se Area GRA SNS (ESN USSR! SSS a (Scale re a Sa Maal ead, WA 
: a = — 
aa ee ee a 
HYDRONIC WITH VARIATION OF T3 PF itons ae a aes) 
pa eg reer ee SAME (PAR SOE LAS ENS REIS Sees GOOSEN SAGEM UE ERE 
See ee arse eee —— 
een eee eee ++ an Pee Sree 
CFM |MDOTAIR|  G _ RE | JH H-AIR CWTR NTU CR | Qa | Tf | 64 
| 750 | 3371 | 1.16 | | 1211 | 00056 | 7.0 | 809 || 4876 em [tooow | Oa [ona aeten | afte 
aia ee a i een aie ey Sec eee Ieee 
COIL DIMENSIONS AIR PROPERTIES Me ot STEAM PROPERTIES ~ ee See Cae 
| COILDIMENSIONS | | —CAXIR PROPERTIES | 
REED ASAIR OCR Ra ae ee ea 
116 | 888 OSE Se See 1.23E-05/ 0.0148 | 0.72 | -30 | INF. i = 
CFM Mee Gt RE TM | UAIR MDOTST -SSTM__GMAK TU | oR ER ee BL 
1.46 | | 1211 | 0.0056 | 7.0 | >o | 658 | 809 | CC] S890 | INF. | 218408 { 0.30 [| 0.00 | 0.26 | 56155 | 39 | 
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DUCT HEATER 2-25- ee Oe ee ia Si: REA A I SR 
anes ASC SS Seat meer OCA Sam Mma nLese Saree RAMEE aA See Se See eae 
or DRAENOR |———| ——|- ——~ anv ropes} ——/-}-——_} — ae See ee eet a a ca 
GOI DIMENSIONS {| AIRPROFERTIES | WATER PROPERTIES Jf Pad 
a a a a aia pe 
2868 | 0312 | | 0.0740 | 0.24 | 4.20605) 0.0148 | 072 | 40 | | 1.002 anes um a a! Ae 
ee ee i ie ee (eared eves RAO ANON IRS: SE 
ceM MDOTAIR) “GRE AR CAR : NTU a a 
i 0.23 | 0.13_| 22526 | 60_| 185 
[4250 [sete 2.06) +(| 4361 | 0.005 | 111 | 4053 | 90 | 1348 aa7e| 4006 | atezwn [018 [ozs {ove {amie {on aan 
[4450 [e516 | 390) —~(| 8047_| b.00a7 | 124 [4083 | 07 | is64_| 4876 | 4886 | 234687 | 0.12 | 0.32 185 
reso | rans | 272 | eras | agoee | as 10.2 4876 | 4886 | 206044 | 0.11 | 0.96 | 0.10 | 27005 | 66 | 184 
P1850 | e314 | 3.05 4053 10.7 | 1995] | ~4876 | 4866 | 200300 | 011 | 0.41 | 0.10 | 2ea16 | 56 | 1e4_| 
2100 9437 | 346 | ~(| 7309 | 0.0030 | 14.5 111 | 2265 456 | 339746 | 0.10 | 0.46 | 0.09 | 30485 | 63 | 164 | 
te ane Coma sao fo | a 
LE At Roe Soa LANE NAS ee oe a 

et a cared peed aan (acca 

ff = 
a noe ai re nse ae NE SRE 
a = ae RO Se eR ee eee aaa — 
oe ee [aks NE a ae ae ee se eee flee ee 
ce (MDOT, GP TH H-WIR CAR” cwiR WU | GR | es | a | m2 | 44 
a7i9_|_i73_| | 3664_| 6.0084 4053 | 8.3 | 1132 gene | tovest as | aas {ats {parva | ee | am 
ee ea ase Se A HT SN eeu (aRaR| 
ia perpen inven eget Ren ae Oa cea eee 
ar ce le Pod ede fe 
oS oe es eed 7 comes MES NER od 
A a Lo a STEAM PROPERTIES Po eT cas 
TAFE | ASAIR | DH “RHO | oP | MU TK ik rcp | 73 | We || A a ea 
Ee I aa eraied cal — 
Ede a ee Se eae aa 
a a a A HSTM | JAR T | Sst 3 GR | EPS] Qs] 2 
/Se64_[ 0.0084 | 101 | >>0 | 65 | 1108 | | 330 | INF. | 226408] 018 | 0.00 | 014 | 31440 | 68 
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po Paes | ee aa eee aS 
ae ee Gea SENT STS Se ee — a Se meet 
[eae Cael, econ 5 ail fete! as 
HYpRON | ale el (eer ee 
ae ane eae ee Saeed eee: cane (ee 
COILDIMENSIONS | | AIRPROPERTIES | | ian Taos, 
AFF | ASAIR| DH | | RHO | cP | MU | K {| PR | T1 es eS ie aaee 
7 | 3079 | 0154 | | 0.0749 | 0.24 | 1.23E-05| 0.0148 | 0.72 | 45 —— ——— 
aS SASOION: NOC SONI) CREAN Re ieeet Coma (PONTOON iiieret 
Se ey ee ee ee = ee naa eae 
| CFM (MDOTAIR, Go| | RE | JH | HAIR | HWIR| UAIR | CAIR CWTR NTU | CR | EP 14 
173 | 082 | | 854 | 0.0062 | 54 | 4053 | 45 | 377 54737 | 0.33 | 0.08 28 | 16152 | 85 | 187 
i910 | 099 | | 1034 [ 0.0058 | 62 | 4053 458 0.31 | 0.09 17313 | __ 83 
| 600 | 2247 | 117 | | 1216 | 0.0056 | 7.0 | 4053 | 56 | 539 78196 | 0.29 | 0141 | 025 | 19158 | 81 
S75 | 2684 | 1.34 | | 1308 | 0.0051 | 7.4 | 4053 | 59 | 620 g9925 | 0.26 | 013 | 0.23 | 20385 | 78 | 186 | 
650 | 2921 | 1.2 7.8 | 4053 | 6.1 701 4886 654 | 0.24 | 014 | 021 | 21277 | 75 | 186 | 
| 700 | 3146 | 1.63 | 8.1 4053 | 6.3 755 4886 teeera 0.28 1018 {020} 22087 | _/_ 188 
Haier eae IEC gney (NR Seen y (eceererey ooee |eee Se) | aaa eee ee ee aa a 
HYDRONIC WITH VARIATION OF T3 pa fe | ae _ as te 
i aa eee ee ere (Ee (ees Te ae | ieee 
es SL SR SE ee Ea a 
fee ae [ose tdl | a es a 
| SFM woOTAR ee JH H-WTR MR NTU | CR | EPs | Q | 72 1 74 | 
ae | 1878 | 882 | ____ft_ 0.0062 | 5.4 | 4063 | 45 | 377 4986 | 60399 | 0.33 | 008 | 0.28 | 16720 | 89 | 202 | 
Lea) pp Heo as) ead ie Se aE aes Eas (epee ee 
i ee ae ee eee ee Eat 
ants sSetere: RE 
AIR PROPERTIES ; STEAM PROPERTIES eee | 
CP MU K | PR | 13 | HFG eee 
| 77 | 3979 | 0.154 | 1.236-05/ 0.0148 | 0.72 INF. 952 po te 
pons ean eee Weel ee 
CFM |MDOTAIR  G | RE | gH | HAIR UAIR CR | EPS | Q | 72 
| 350 | 1573 | 0.82 5.4 47 INF. | 73612 | 0.34 [0.00 | 0.29 | 21428 
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APPENDIX E. MANUFACTURER’S DATA FOR 
LARGER CAPACITY FANS COMPATIBLE WITH 
B25 UNIT HEATERS 
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The 


New York Blower 


Company® 





7660 QUINCY STREET—WILLOWBROOK, ILLINOIS 60521 






: - be} n st Be 7 e vi fe 2 
¢ : 4 ‘ 3 e of | \, os = 


Os tes Pas sity yt st . 
Model Max* 


SN82-H 20 1550 442 316 270 

SN82-H-3_ | t¥20 | 1550/1300/1100 | 442/371/314 } 316/265/224 | 270/226/192 

SN102-H 20 1550 870 755 720 
SN102-H-3} f¥20 | 1550/1300/1100 | 870/730/617 | 755/633/536 | 720/604/511 


1150 920 850 = 




















































1475 
1475/975 


SN122-H 
SN122-MH 


1815 
1815/1200 


1675 
1675/1106 


1650 
1650/1090 





1725/1140 










































SN142-M V2 1350 1160 1100 _— — 
SN142-H 2100 1990 1960 1840 1680 
SN142-MH 1725/1140 2100/1390 1990/1315 1960/1295 1840/1216 1680/1110 


































Vg 
V4 
V4 
Yo 
iy 
V4 
lo 

v4 


SN162-M 1140 2000 1800 1750 1450 
SN162-H 1750 2950 2830 2800 2650 2500 
SN162-MH 1725/1140 2900/1915 | 2790/1840 2760/1756 2600/1718 2440/1610 






































SN182-M a 1140 2610 2400 2340 1960 _ — 
SN182-H Vi 1725 3920 3750 3700 3490 3280 3000 
SN182-MH| ¥ 1725/1140 3920/2590 3750/2480 3700/2440 3490/2305 3280/2165 3000/1980 





SN202-M V4 1140 3570 3260 3200 - 2810 — — —_ 
SN202-H ¥4 1725 5300 5100 5000 4820 4600 4350 92 
SN242-L VY 1140 4400 4150 4080 3700 — — — 
SN242-M Y2 1140 5380 5100 5030 4650 4200 _ 2 
SN242-H | *3/4 6020 5600 5120 4480 9 


* Maximum BHP over cataloged range. Motors are rated on internal temperature rise rather than 
NOTE Static pressure rating on nameplate HP. 


pultiapeed tens i> atthe higher t Shaded-pole motor. Three-speed capacities shown are obtainable with 3-speed switch furnished 
with unit. 


* Available in 3-phase only. 














speed. Low speed capacities are 
shown for the indentical system. 
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DIRECT-DRIVE PROPELLER FANS 
EXHAUST OR SUPPLY 


® Eight wheel diameters—8" through 24". 

@ 243 to 6420 CFM— up to 1/2" static pressure. 

® Panels—square stee! construction with streamlined venturi 
inlet...venturi is reversed in supply-fan panels. ..baked-green 
alkyd finish. 

@® Wheels—aluminum blades with steel hubs. 

® Motor mounts—wire-guard-type motor mount (see photo at left] 
is standard on all direct-drive units...guard is zinc-plated steel. 

@ Motors—standard motors are totally enclosed air over with pre- 
lubricated ball bearings except 12 and 20 HP motors, which 
are shaded-pole totally enclosed permanently lubricated sleeve- 
bearing type. 
Motors 4 HP and larger are suitable for either horizontal or 
vertical service...specify “for vertical mounting” to have wheel 
locked to motor shaft...1/20 and 1412 HP motors are not suitable 
for vertical service. 





b}-—H-—--| 
H = min. wall 
thickness 





Mounting hole . 
“ag 


tbs. 
8-6 25 









shutter 








not ae 
available 
80 
95 
t E and G based on longest motor used for each size fan. * Shipping weights shown are maximum and include totally enctosed motors and weight of packaging. 
NOTE: Exhaust units are avaitable with either automatic or motorized shutters. Supply units require motorized supply shutter. 
When ordering, specify compfete mode! number as shown on page 3. Dimensions not to be used for construction unless certified. Tolerance: + Ve" 
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APPENDIX F. COMPUTATION OF AIR SIDE 
CONVECTION COEFFICIENT USING 
MANUFACTURER’S DATA 
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Flom 400 a ; 
=> Vai = 750 er a Set = 5 00 ‘Ye i 
Ag, ~ ( lx]. 5 )et* 


m4 (2 [2 evtsfnch | nod SOO lGore 
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Q= G, AT = mee, aT = CAVG aT 
c f “7 nar QT. \¢ GOmIn 
= lo.o7t4 pi WE) he) (oe A) 712€) ST 


= G2449 8h, 








Ee ee tes * tee 
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“ere =e eon 


=> NTV= (1-2) 
: Jn (12 034) 
- 0.42 
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